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The  3D-BSW  Model  Applied  to  Climatology 

of  Small  Areas  and  Lines 


1.  INTRODUCTION 

A  previous  report1  describes  the  procedures  for  the  simulation  of  a  spatial 
set  of  conditions  by  use  of  the  three-dimensional  Boehm  Sawtooth  Wave  Model 
(3D-BSW).  The  present  report  describes  the  use  of  the  model  for  estimating  the 
climatic  probability  of  spatial  events  by  algorithms  based  on  repetitive  simulations. 

Special  attention  is  given  to  cloud-cover  probabilities.  Fields  of  clouds  are 
simulated  stochastically,  incorporating  areal  sizes  and  lineal  separations  of  the 
clouds,  both  in  the  horizontal  and  in  the  vertical.  Using  these  simulations, 
synthetic  climatology  is  developed  to  estimate  the  likelihood  of  events  such  as 
an  overcast  of  clouds  covering  a  floor  space  of  varying  dimensions  or  the 
probability  of  a  clearing  in  the  clouds.  One  special  problem  addressed  is 
the  probability  of  a  cloud-free  interval  along  a  line  of  travel. 


(Received  for  publication  7  August  1987) 

1.  Gringorten,  1. 1.  (1984)  A  Simulation  of  Weather  in  3D  Space. 
AFGL-TR-84-0267,  ADA  155221. 
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2  3  4 

The  approach,  as  in  previous  reports,  ’  '  is  to  develop  a  model  that  will 
depict  a  synoptic  situation  that  statistically  resembles  the  natural  state  of  the 
weather.  In  many  repetitive  applications  of  the  model  the  correlation  of  events 
at  points  of  known  distance  apart  should  resemble  the  natural-state  correlations. 
Additionally,  a  multitude  of  applications  of  the  model  should  yield  frequencies  of 
events  that  will  approximate  their  climatic  probabilities. 

In  the  previous  report  a  set  of  isopleths  of  rainfall,  in  one  actual  24 -hr 
period  in  New  England,  was  shown  as  an  example  of  the  kind  of  field,  with  all  its 
variability,  that  could  be  simulated  by  a  model.  The  3D-BSW  model  is  the  basis  for 
a  Monte  Carlo  method  to  produce  climatological  probability  estimates.  To  illus¬ 
trate  this  purpose.  Figure  1  gives  such  estimates  of  the  probability  of  24-hr 
rainfall  in  New  England,  in  reasonable  agreement  with  the  results  of  surveying 
January  24-hr  amounts,  in  10  years  at  some  140  stations  through  New  England. 
As  Figure  1  shows,  the  single-station  frequency  of  no  rain  is  63  percent;  that  of 
rain  not  exceeding  13  mm  U/2  inch)  throughout  Massachusetts  is  84  percent,  and 
is  accurately  estimated  as  such  by  the  model.  The  3D-BSW  model  extends  the 
estimates  to  varying  areal  sizes,  and  shows,  for  example,  that  in  1000  km  the 

probability  of  no  rain  is  57  percent.  Or,  the  probability  of  rain,  of  any  amount, 

2  2 
somewhere  in  1000  km  is  43  percent.  Also,  somewhere  in  the  1000  km  area 

there  is  an  8  percent  probability  of  a  13-mm  (1/2  inch)  rainfall,  and  so  on. 

(One  exception  to  the  application  of  the  model  is  on  Mount  Washington,  N.  H. , 
where  the  precipitation  is  heavier  than  in  the  rest  of  the  region.  ) 

As  in  preceding  reports,  model  development  is  done  with  the  equivalent  normal 
deviate  (END),  symbolized  as  y,  of  the  meteorological  variable  (X).  The  trans¬ 
formation,  back  and  forth,  between  X  and  y,  is  accomplished  through  the  cumulative 
probability,  Pr(X)  or  Pr(y).  Symbolically, 

X* — ♦  Pr(X)  =  Pr(y)* — *y  .  (1) 

Several  methods  of  accomplishing  this  transformation  have  been  mentioned  and 
illustrated  in  previous  reports.  Whatever  method  is  adopted,  a  one-to-one  relation 
will  exist  between  y  and  X,  or  between  y  and  Pr(X). 


2.  Gringorten,  1.  1.  (1979)  Probability  models  of  weather  conditions  occupying  a 

Line  or  an  area,  J.  Appl.  Meteorol. ,  18:957-977. 

3.  Whiton,  K.C.,  Berecek,  E.  M. ,  and  Sladen,  J.G.  (1981)  Cloud  Forecast 

Simulation  Model,  USAFETAC  Scott  AFB,  IL  62225,  USAFETAC/TN-81-004. 

^126  pp. 

4.  Burger,  C.F.,  and  G  ringorten,  1.  1.  (1984)  Two-Dimensional  Modeling  for 

Lineal  and  Areal  Probabilities  of  Weather  Conditions,  AFGL-TR-84-0126, 
ADA  147970,  58  pp. 
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Figure  1.  Showing  the  Frequency,  at  the  X's  of  24-hr  January 
Rainfall  at  a  Typical  New  England  Station  (on  the  left-hand  ^ 
ordinate);  of  the  Maximum  in  Massachusetts,  Area  21,350  km  ; 
and  in  the  Six-State  New  England  Region,  Area  169,  824  km  . 
The  solid  curves  are  found  by  the  application  of  the  3D-BSW 
model,  with  scale  distance  of  6.  22  km.  They  show  the 
probability  of  no  rain,  of  rain  less  than  13  mm  (1/2  inch)  and 
rain  less  than  50  mm  (2  inches),  throughout  an  area  (A) 
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Throughout  this  work  we  assume  the  weather  is  stationary  stochastic.  While  there 
is  much  variation  in  the  condition  (X)  from  place  to  place,  with  level  of  elevation 
and  in  time  of  day  or  season,  it  is  assumed  that  the  "true"  average,  median  or 
mode,  and  the  true  probability  distribution  of  values  remain  unchanged  with  time 
or  era.  Although  this  premise  conflicts  with  the  occasional  declaration  of  changing 
climate,  it  is  acceptable  in  a  human  lifetime  of  climatic  events,  including  extremes. 

The  problems  and  solutions,  as  presented  herein,  might  also  be  viewed  as  a 
generalization  of  the  often-asked  questions  about  joint  frequencies  or  probabilities 
of  events,  such  as  simultaneous  cloud  cover  at  two  or  more  stations.  This  study  deals 
with  the  aggregate  of  events  throughout  a  specified  area,  answering  such  questions  as 
the  likelihood  of  consistently  high  values  throughout  the  area  or  in  a  measurable  fraction 
of  the  area.  That  is.  the  3D-BSW  model  is  extended  to  estimate  the  probability  of 
events  occurring  in  an  area  or  given  fraction  of  an  area,  or  along  a  line. 


2.  THE  3D-BSW  MODEL 

With  the  3D-BSW  model,  a  value  of  the  END  (y)  is  obtainable  at  every  point 
in  3D-space.  Care  must  be  taken  to  modify  the  vertical  measure  (w1)  by 
multiplying  it  by  a  factor  q  (=  50,  say)  to  make  the  modified  vertical  measure  (w) 
comparable  to  the  horizontal  measure  in  the  (u,  v)-plane.  Thus: 

w  =  q.  w*  .  (2) 

2.1  Development 

For  the  3D-BSW  model  each  formation  of  uniform  sawtooth  waves  fills  all  of 
three-dimensional  space.  (The  waves  do  not  move. )  Each  wave  formation  has  a 
constant  wavelength  (A)  that  becomes  a  characteristic  parameter  of  the  model.  The 
wavelength  may  vary  from  tens  of  kilometers  to  hundreds  of  kilometers. 

A  parameter,  called  the  scale  distance,  symbolized  as  r,  is  proportional  to  A. 
In  earlier  work  the  scale  distance  was  defined  as  the  distance  over  which  the 
correlation  coefficient  is  0.  99.  In  the  study  of  clouds,  r  was  found  to  be  approxi¬ 
mately  1  km  when  A  was  several  hundred  kilometers.  In  the  3D-BSW  model,  we 
set  the  ratio  of  A  to  r  at 

A  =  256  r  km  .  (3) 

This  equivalency  relationship  permits  us  to  use  the  scale  distance  of  the  earlier 
2 

Gringorten  Model  B  as  the  scale  distance  of  the  3D-BSW  model  (see  Section  2.4.  1 
for  amplification). 
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2.  1.  1  ELEMENTS  NEEDED  FOR  A  SINGLE 
WAVE  FORMATION 

Consider  one  wave  formation  in  the  3D-space,  oriented  so  that  the  flat  surface 
or  plane  that  contains  one  constant  phase  is  parallel  to  the  plane  surface  ABC 
(Figure  2).  Let  X(u,  v,  w)  be  a  point  on  the  plane  surface  (ABC)  with  coordinates 
(u,  v,  w).  Let  P  be  the  closest  point  on  this  plane  to  the  origin  0(0,  0,  0).  Then  the 
line  OP  will  be  perpendicular  to  the  plane,  oriented  by  angles  (a.  ft  ,y  )  between 
OP  and  the  three  axes.  If  the  length  of  OP  is  D,  then  from  solid  geometry  the 
equation  of  the  plane  containing  P  and  X  is 


u.  cos  a  +  v.  cos  /3  +  w,  cos  y  =  D  . 


(4) 


Let  X  be  the  angular  measure  from  the  U-W  plane  to  the  vertical  plane  containing 
OP.  (If  the  W-axis  were  the  earth's  axis,  X  would  be  longitude.  )  From  spherical 
triangles  (Figure  3), 


cos  a  =  sin  >  .  cos  X 


cos  /3  =  sin  y  .  sin  X 


from  which 


(5) 


2  2  2 
cos  a  +  cos  0  +  cos  y  =  1 


(6) 


Or,  one  of  (a,  j3 ,  y  )  is  determined  when  the  other  two  are  given. 

Let  the  leading  edge  of  the  middle  wave  of  a  wave  formation  miss  the  origin 
by  a  distance  H,  and  define 


h  =  H/A  . 


(7) 


There  are  three  degrees  of  freedom  in  the  random  choice  of  a  wave  formation. 
We  choose  h  by  random  process,  thus  using  one  of  the  three  degrees  of  freedom: 


OShSl.O  , 


(8) 


h  being  uniform  on  (0,  1). 

The  remaining  two  degrees  of  freedom  must  be  used  to  obtain  values  for 
(at,  &  ,y  ).  We  have  two  alternative  procedures  for  doing  this.  While  one  is  a  more 
direct  procedure,  the  other  will  save  computer  time. 


Figure  2.  Illustration  of  a  Plane  (ABC).  The  point  (X)  is  considered 
to  be  on  this  plane  surface.  Point  (P)  is  such  that  OP  is  the  shortest 
distance  from  the  origin  to  the  plane,  and  has  angles  (a,  /3  ,y  )  with  the 
U-,  V-,  W-axes.  The  projection  of  OP  onto  the  U-V  plane  has  angle  \ 
with  the  U-axis 


2.  1.  2  THE  LAMBDA  METHOD  FOR  OBTAINING 

(a,  1 3  ,  y  ) 

Consider  the  sphere  with  radius  (OP).  (Figure  3  shows  one  octant.)  The  shell 
2 

has  a  total  area  of  4jtD  .  An  element  of  the  area  (dA)  is  given  by: 


Figure  3.  Illustration  erf  One  Octant  of  the  Sphere  Centered 
at  the  origin  (0).  The  line  (OP)  of  length  (D)  has  angular 
sides  (a.P  ,y)  with  the  three  axes  (U,  V,  W)  and  lies  in  the 
plane  that  makes  angle  (X)  with  the  (U-W)  plane 


dA  =  (D.dy).  (D.  siny.dX) 

(S 

=  D2  .  d (-cos  y  )  .  dX  . 

In  the  random  selection  of  the  wave  formation,  the  point  P  should  have  the  same 
probability  of  appearing  in  one  element  of  area  (dA)  as  in  any  other  on  the  sphere 
Thus  we  choose  two  random  numbers  (R^,  R^),  each  between  0  and  1,  and  set 


cos  y  *  2R  - 1,  to  make  -1.0  i  cos  y  i  1.0 
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where  cos  y  is  uniform  on  (-1,  +1) 

A  =  2jtR^,  to  make  0.0  5  A  %  2ir  . 
where  A  is  uniform  on  (0,  2n) 

Equations  (5)  give  cos  a  and  cos  0,  which,  together  with  cos  y  define  D 
[Eq.  (4)]  for  each  point  (u,  v,  w)  in  the  3D  space.  We  will  refer  to  this  as 
Alternative  (1). 

2,  1.3  THE  ACCEPTANCE -REJECTION  METHOD  FOR 
OBTAINING  (a,  0  ,  y  ) 

Choose  three  random  numbers  (R^,  Ry>  Rw>  between  0.0  and  1.0.  and  set 

r  =  R  -  1/2  to  make  -1/2  s  r  §  1/2  . 
u  u  u 

ry  =  Ry  -  1/2  to  make  -1/2  =  r  S  1/2,  (11) 

rw  =  Rw  -  1/2  to  make  -1/2  =  rw=  1/2  . 

The  point  with  coordinates  (ru>  r^,  r^)  will  lie  within  a  unit  cube.  To  assure  that 
any  one  direction  of  the  line  from  the  origin  to  (r^,  ry,  r  )  will  be  as  likely  as  any 
other  direction,  we  impose  the  restriction: 

(OP')  =  Jr2  +  r2  +  r2  s  1/2  .  (12) 

U  V  w 

If  the  length  (OP1)  is  greater  than  (1/2),  then  we  reject  the  three  random  numbers 
and  try  another  group  of  three  random  numbers.  The  probability  that  the  set 
(Ru«  Ry,  Rw>  will  be  rejected  will  be  (1  -  4 /3w(  1  / 2)^]  =  0.476.  On  the  average, 
this  alternative  procedure  requires  the  generation  of  (3/0.  524)  =  5.  7  random 
numbers  per  wave  formation  for  orientation  plus  an  additional  random  number 
for  phase. 

With  the  first  acceptable  set  of  random  numbers,  we  find 
cos  a  *  r^/ (OP')  , 

cos  0  =  rv/(OP’)  ,  (13) 

cos  y  =  rff/(OP')  , 
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which  then  define  D  [Eq.  (4))  for  each  point  (u,  v,  w)  in  the  space.  We  will  refer 
to  this  method  as  Alternative  (2). 

2. 1.4  COMPARISON  OF  ALTERNATIVES  FOR 

GENERATING  WAVES 

While  Alternative  (2)  requires,  on  the  average,  the  generation  of  6.7  random 
numbers,  as  opposed  to  three  for  Alternative  (1),  it  does  not  require  the  deter¬ 
mination  of  a  sinusoidal  function,  that  consumes  more  computer  time  than  finding 
the  random  numbers.  There  is  another  attractive  feature  in  Alternative  (2):  it 
can  be  extended  to  higher  dimensions  than  the  three  for  space,  to  include  an  addi¬ 
tional  dimension  for  time  or  to  introduce  additional  influences  on  the  wave  forma¬ 
tion.  (This  will  be  explored  in  later  reports. )  However,  care  must  be  taken  to 
eliminate  the  selection  of  absolute  values  for  (r^,  r^,  rw>  that  will  make  (OP')  too 
small  for  computer  operation. 

2.1.5  WAVE  HEIGHT 

The  wave  height  (x)  of  one  wave  formation,  at  a  point  (u,  v,  w)  is  given  by 

x  =  FRA  (h  +  D/A)  (14) 

where  FRA  denotes  the  fractional  part  of  the  bracketed  quantity.  While  x(u,  v,  w) 
is  termed  "wave  height",  it  is  difficult  to  visualize  in  3D  space.  For  our  purpose, 
it  is  acceptable  as  a  property  of  the  wave  formation  at  the  point  (u,  v,  w),  and  is 
defined,  mathematically  and  unequivocally,  by  Eq.  (14).  Where  x'  =  h  +  D/A,  x' 
can  be  negative.  The  symbol  INT(x')  is  used  for  the  integer  part  of  x'.  If,  as  in 
FORTRAN,  this  is  taken  to  be  the  integer  part  with  its  sign,  then  the  wave  height 
(x)  at  a  point  (u,  v,  w)  is  given  by 

x  =  x'  -  INT(x')  when  x'  i  0  , 

=  x'  -  INT(x')+  1  when  x'  <  0  . 

If,  as  in  BASIC,  the  symbol  INT(x')  is  taken  to  be  the  next  lesser  whole  number, 
then,  simply,  the  wave  height  is  given  by 

x=  x'  -  INT(x')  for  all  x'  . 

The  wave  height  (x)  varies  uniformly  from  0  to  1.  0  in  the  3D-BSW  model. 
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2.  1.6  GENERATING  MULTIPLE  WAVE  FORMATIONS 

The  above  procedure  produces  one  wave  formation  and  one  set  of  values  for 
cos -a,  cos  (i ,  cos  y ,  and  h.  Repeated  L  times,  there  will  be  L  wave  formations, 
with  L  wave  heights  (X|,  l  =  l(l)L)  at  (u,  v,  w)  whose  average  wave  height 

L 

x  -  £  x.  (u,  v,  w)  (15) 

L  i=l  1 

will  have  a  mean  of  (1/2)  and  variance  of  (1/12L),  and  will  have  an  approximate 
normal  (Gaussian)  distribution.  For  sufficiently  large  L,  the  END  (y)  at 
X(u,  v,  w)  is  given,  practically,  by 

y  (u,  v,  w)  =  nTTSTT.  (x  -  1/2)  .  (16) 

The  END  (y)  can  thus  be  calculated  for  any  point  (u,  v,  w)  in  the  3D  space. 

If,  for  convenience,  L  is  made  equal  to  12,  then 
12 

y  (u,  v,  w)  =  Zj  x.  (u,  v,  w)  -  6  .  (17) 

i=  1  1 

For  a  single  simulated  field,  a  minimum  of  36  random  numbers  using  Alternative  (1), 
or  by  Alternative  (2)  an  expected  maximum  of  80  random  numbers,  are  needed  to 
generate  all  values  of  the  END  (y)  at  all  points  in  the  three-dimensional  field. 

22,  Simulation  or  Depiction 

To  obtain  a  picture  in  a  horizontal  surface,  within,  say  a  square  200  km  to  the 
side,  we  keep  w  constant  (conveniently  zero),  and  compute  y(u,  v,  w  =  0)  for  as 
many  points  (u,  v,  w=  0)  in  the  square  as  we  wish.  Likewise  we  can  generate  a 
vertical  cross-section  by  keeping  u  constant,  or  v  constant,  or  by  keeping  a  linear 
combination  of  (u,  v)  constant.  For  an  overall  procedure  to  produce  a  field  of 
values,  see  Appendix  A. 

Samples  of  simulated  cloud  structure  are  given  in  the  previous  report1 
(October  1984),  one  in  a  horizontal  plane  representing  a  cloud  cover,  and  four 
other  examples  of  clouds  in  a  vertical  cross-section. 
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23  Probabilities  -  Description  of  the 
Monte  Carlo  Method 

The  above -described  procedure  produces  END-values  in  all  three  dimensions, 
or  in  the  two  dimensions  of  a  horizontal  surface  or  vertical  cross-section,  or  along 
a  line  of  travel  or  line  of  sight.  Suppose  a  field  of  END-values  (y)  is  generated  in 
a  horizontal  square  of  area  A.  Among  the  many  values  of  y  in  the  square  there  will 
be  a  maximum,  yma3£»  also  a  value  that  will  be  the  (9/10)ths  maximum,  also 
(8/10)ths  maximum,  and  so  on.  Taking  any  single  value,  yQ  (say  =  0.4),  as  a 
threshold,  there  will  be  a  certain  probability  that  it  will  be  the  maximum  in  the 
area  and  will  not  be  exceeded  in  the  entire  square.  There  will  be  a  higher  prob¬ 
ability  that  it  will  not  be  exceeded  in  all  except  (l/10)th  of  the  area,  still  higher 
in  all  except  (2/10)ths  of  the  area,  and  so  on.  Lastly  there  is  the  probability 
that  is  1.0  minus  the  probability  that  it  will  be  exceeded  everywhere,  when  it  will 
be  the  minimum  in  the  square.  The  areal  climatology  of  a  region  should  consist 
of  information  on  such  probabilities.  The  same  kind  of  climatology  could  be 
obtained  for  events  or  conditions  along  a  line. 

At  this  stage  we  must  admit  to  an  unfulfilled  wish.  An  analytical  solution  is 

desirable  to  determine  the  probability  distribution  of  the  maximum  in  a  space, 

similar  to  the  solution  for  the  maximum  in  an  Ornstein-Uhlenbeck  temporal 
5  6 

process.  *  Not  having  such  an  analytical  approach,  we  have  resorted  to  approxi¬ 
mate  solutions  by  the  Monte  Carlo  method. 

2.  3. 1  MULTIPLE  FIELDS  OF  ENDS 

In  an  area  (A)  or  on  a  line  (L)  a  finite  number  of  END-values  can  be  found  by 
the  above-described  procedures.  For  a  single  such  field,  a  maximum  value  of  y, 
the  above  percentiles  and  minimum  value  can  be  determined.  The  greater  the 
number  of  y-values  that  are  calculated,  the  closer  will  be  the  maximum  to  the  true 
maximum  of  the  area  or  the  line.  Now,  let  us  generate  many  such  fields  or 
"snapshot"  pictures  (M  =  25,000,  say)  of  area  (A)  or  line  (L).  There  will  be  M 
trial  values  of  the  maximum  in  A  or  L.  These  values  collectively  constitute  a 
frequency  distribution  of  the  maximum,  and  may  be  regarded  as  a  good  approxi¬ 
mation  to  the  true  probability  distribution  of  the  maximum.  The  same  will  be  true 
for  all  percentiles  and  for  the  minimum. 


5.  Keilson,  J. ,  and  Ross,  H. F.  (1975)  Passage  time  distributions  for  Gaussian 

Markov  (Ornstein-Uhlenbeck)  statistical  process.  Selected  Tables  in 

Mathematical  Statistics,  3:233-237. 

"  "  •  ■  ■  '  '  «■» 

6.  Gringorten,  I.I.  (1982)  The  Keilson-Ross  Procedure  for  Estimating  Climatic 

Probabilities  of  Duration  of  Weather  Conditions,  Afgl-T^-82-(Tii6. 

ADA  119860.  - - 
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2.3.2  TYPES  OF  PROBABILITIES  DEFINED 


For  clarification  in  the  following  development,  the  symbol  P^  is  used  for 
probability  of  exceedance  of  the  value  (X)  at  a  single  point.  The  symbol  PQ,  as 
used  herein,  is  the  cumulative  probability,  and  is  the  probability  that  X  will  not 
be  exceeded  at  a  single  point.  That  is, 

P  =  (1  -  P')  .  (18) 

o  o 

It  is  assumed  that  PQ  is  the  same  for  any  point  in  the  area  or  the  line.  The 
climatic  probability  of  the  weather  condition  defines  Pq,  which,  in  turn,  defines 
its  corresponding  END  (yQ).  As  mentioned  in  Section  1  there  is  a  one-to-one 
correspondence  between  the  END  (y)  and  the  meteorological  variable  (x).  A 
threshold  (X)  such  that  x§  X  has  its  corresponding  value,  y  ,  such  that  y  =  yQ. 

The  single-point  END  (yQ)  will  be  negative  when  Pq  <  1/2,  positive  when 

P^  >  1/2.  (For  sky  cover  Pq  is  the  cloud-free  probability  in  a  single  direction.) 

The  symbol  PA  will  be  used  for  the  probability  of  a  fraction  or  less  coverage  of 

2 

an  area.  An  example  is  the  probability  of  3/10  or  less  cloud  cover  in  1000  km  . 
Coverage  can  also  apply  to  the  threshold  of  a  continuous  variable  such  as  tempera¬ 
ture.  An  example  is  the  probability  that  0.  25  of  the  area  has  temperatures  below 
freezing. 

The  symbol  PL  will  be  used  for  the  probability  of  a  fraction  or  less  coverage 
of  a  line  segment,  for  example,  the  probability  that  along  the  route  between 
New  York  City  and  Washington,  D.  C.,  rain  would  occur  0.33  or  less  of  the  distance. 
Coverage  is  without  regard  to  whether  the  rain  occurred  in  one  unbroken  interval 
or  many  small  intervals  as  long  as  their  sum  is  0.  33  or  less  of  the  total  distance. 

The  symbol  PI  will  be  used  for  the  probability  that  at  least  one  interval  some¬ 
where  along  a  line  segment  has  weather  continuously  below  the  threshold  (yQ). 

For  example,  suppose  a  helicopter  flying  through  patchy  fog  can  visually  navigate 
if  no  fog  patch  is  greater  than  1.5  km.  What  is  the  probability  that  there  are  one 
or  more  fog  patches  greater  than  1.5  km  anywhere  along  a  200-km  route?  A 
weather  condition  can  be  a  specific  kind  of  weather  (for  example,  thunderstorm)  or 
all  values  exceeding  (above  or  below)  a  specific  threshold  (for  example,  pressures 
below  992  mb).  Note  that,  when  the  interval  is  as  long  as  the  total  line  segment, 
then  PI  is  the  same  as  PL  for  zero  coverage. 

A  line  interval  (s')  may  be  visualized  as  a  window  sliding  along  a  longer  line 
of  travel  (T);  at  any  position  of  the  window,  within  s'  there  will  be  a  maximum  value; 
the  lowest  of  these  maximum  values  is  the  threshold  condition  whose  probability 
(PI)  we  wish  to  find. 


12 


2  A  Areal  Coverage  (PA) 

The  goal  of  this  section  is  to  find  a  satisfactory  approximation  for 
PA  (X,  A,  F,  r)  =  PA(yo,  A,  F,  r) 

which  is  the  probability  that  a  weather  condition  (x)  will  not  exceed  the  threshold 
(X)  (or  the  END  (y)  will  not  exceed  yQ),  except  in(F/10)ths  orless  of  the  area  (A), 
given  the  scale  parameter  (r).  That  is,  1-PA  (yQ,  A,  F,  r)  is  the  probability  that 
yQ  is  the  (F/10)  highest  value  of  y  in  the  area  (A). 

If  F  =  0,  then  the  threshold  weather  condition  is  exceeded  nowhere  in  the  area. 
For  example,  what  is  the  probability  of  no  clouds  anywhere  in  New  Jersey,  that  is, 
the  state  is  completely  clear? 

If  F=  10,  then  the  threshold  weather  condition  is  exceeded  everywhere  in  the 
area.  For  example,  what  is  the  probability  the  relative  humidity  is  above 
80  percent  everywhere  in  Iowa? 

If  0  <  F  <  10,  then  PA  is  the  probability  that  the  threshold  weatner  condition 
is  not  exceeded  except  in  a  fraction  of  the  area  less  than  or  equal  to  (F / 1 0). 

2.4.1  STANDARDIZED  UNITS 

At  this  stage  of  the  development,  we  need  to  standardize  our  units.  A  square, 
2 

of  area  A  (km  ),  has  side  s'  where 


s'  =  nTa  km 


(19) 


The  standardized  distance  (s)  is 

s  =  s'/r  (20) 

where  r  is  the  scale  distance.  Or 
r  =  s'/s  km  . 

The  scale  distance  (r)  had  been  defined,  in  earlier  work,  as  the  distance  over 
which  the  correlation  coefficient  is  0.  99.  In  the  study  of  clouds,  (r)  was  approxi¬ 
mately  1  km  when  the  wavelength  was  found  to  be  several  hundred  kilometers. 
Arbitrarily  we  have  chosen,  for  standardization,  the  ratio  of  A  to  r  at  256:1.  That 
is,  A  =  256  r  km  in  the  3D-BSW  model,  making  the  two  parameters  (A  and  r) 
alternatives  to  each  other. 
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In  previous  work  with  New  England  24-hr  precipitation,  we  found  the  scale 
distance  equal  to  5  to  10  km  (6.  22  km  in  the  example  of  Figure  1).  For  radar 
echoes  on  a  PPI-scope  the  scale  distance  ranged  from  1  to  2  km  in  summer, 

2-1/2  to  4-1/2  km  in  winter.  The  following  sections  will  present  typical  results 
for  sky-cover  data. 

For  convenience,  the  term  (z)  is  introduced  to  serve  as  a  dimensionless 
substitute  for  A  or  s',  such  that 

z  =fn(s'/r)/fn2=fnWA/r)/<n2  (21) 

or 

s  =  2Z  where  s  =  s'/r  (22) 

or 


With  the  parameter  A  or  r  given,  there  will  be  a  one-to-one  correspondence 
between  z  and  A.  Symbolically 

PA(yo,  A,  F,  r)  =  PA  (yQ,  z.  F)  . 

2.4.2  GENERATING  MAPS  OF  END-VALUES 

A  computer  program  was  compiled  to  provide  the  Monte  Carlo  estimates  of 
probabilities:  PA  (yQ,  z,  F)  for  yQ  =  -4.  0(0.  5)4.  0,  z  =  0(1)8,  and  F  =  0(1)5. 
Corresponding  to  the  nine  choices  for  z,  nine  square  areas  were  selected  for  the 
study,  that  had  standardized  sides: 

s  =  2Z  =  1,  2.  ...  256 
for  which 

A  =  1,  4,  .. ..  (256)2  for  r=  1)  . 
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By  the  method  described  on  the  previous  page,  and  as  described  in  the  previous 

report,  a  field  of  y-values  was  generated  to  simulate  each  map  of  a  total  of 

25,000  maps.  (In  applying  Eq.  (14)  the  parameter  was  set  at  A  =  256  km.)  Each 

2 

square  was  gridded  to  yield  (21)  =441  grid  points  at  which  the  y-values  were  found, 
which  were  then  ordered  by  size.  The  largest  value  was  assigned  to  the  distribu¬ 
tion  for  y(max).  The  (np)th  largest  value  from  the  top,  where 

np  =  TRUNC  }  44 IF /  10  }  +  1,  F  =  1(1)5  (23) 

was  assigned  to  the  distribution  of  the  (F/10)th  highest  value  of  y  in  the  area  (A). 

The  (F /10)th  lowest  y-value  should  be  comparable  to  the  (F/l0)th  highest  with 
the  sign  of  y  reversed.  Another  y-value,  therefore,  was  obtained  by  selecting  the 
(n^'lth  value  of  y  where 

np'  =  441  -  TRUNC  |441F/10|  ,  F=4(-l)0  (24) 

and,  after  changing  the  sign  of  y,  assigning  it  to  the  distribution  of  the  (F/10)th 
highest  of  y  in  the  area  (A). 

The  six  probability  distributions  of  PA  for  the  six  values  of  F  =  0(1)5  were 
approximated  by  the  frequency  distributions  of  the  y-values  in  the  25,  000  maps. 

The  results,  in  nine  columns  for  the  nine  z's,  in  order  of  increasing  END  (yQ). 
are  shown  in  Tables  1(0)  to  1(5),  each  table  for  one  value  of  F. 

These  values  were  plotted  in  six  charts  [Figures  4(0(  to  4(5)],  and  curves 
were  drawn,  with  smoothing  at  the  extremes.  Each  chart,  one  for  each  value  of 
F  =  0(1)5,  has  probability  along  the  ordinate,  z  along  the  abscissa.  Each  curve 
on  each  chart  is  drawn  for  one  value  of  yQ  [=  -3.  5(0.  5)4.  0].  Each  chart  shows 
the  cumulative  probability  of  yQ  that  will  be  exceeded  in  only  (F/10)ths  of  the  area. 
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Table  1(0).  Estimates  of  the  Probability  Distribution  of  PA(y  ,  z,  F=  0)  Where 

z  2  ° 

yQ  is  the  Maximum  in  the  Area  A  =  (r.  2  )  ,  z=  0(1)8,  r  is  the  Given  Scale 
Distance  (km).  The  sample  size  was  25,  000  maps 


Z 

^o 

0 

1 

2 

3 

4 

5 

6 

7 

8 

-3.5 

0.00009 

0. 00015 

0.00011 

0.00001 

- 

- 

- 

- 

- 

-3.  0 

0. 00097 

0.00101 

0.00081 

0.  0005  9 

0. 00027 

- 

- 

- 

- 

-2.5 

0.00489 

0. 00481 

0.00401 

0.00251 

0.00123 

0.00013 

- 

- 

- 

-2.0 

0.0216 

0.0198 

0.0162 

0.0106 

0.00525 

0. 0075 

- 

- 

- 

-1.5 

0.0611 

0.  0581 

0.0486 

0.0376 

0.  0212 

0.00623 

0.00033 

- 

- 

-1.0 

0.  151 

0.  143 

0.  127 

0.  102 

0.  0643 

0. 0263 

0.  00289 

- 

- 

-0.  5 

0.298 

0.  285 

0.  260 

0.  220 

0.  155 

0.0798 

0.0171 

0.00009 

- 

0 

0.483 

0.471 

0.443 

0.  390 

0.  306 

0.  190 

0. 0637 

0. 00223 

- 

0.5 

0.674 

0.663 

0.637 

0.590 

0.503 

0.  368 

0.  183 

0.0224 

0. 00001 

1.0 

0.828 

0.819 

0.799 

0.760 

0.693 

0.571 

0.374 

0.  110 

0.00157 

1.5 

0. 9298 

0.9234 

0.9132 

0.890 

0.848 

0.765 

0.605 

0.313 

0. 0332 

2.0 

0.9749 

0.9718 

0.9677 

0.9591 

0.9387 

0.897 

0.802 

0.590 

0.227 

2.5 

0.99377 

0.99275 

0.99059 

0.9880 

0.9814 

0.9666 

0.9260 

0.827 

0.598 

3.0 

0.99871 

0.99833 

0.99771 

0.99697 

0.  99511 

0. 99095 

0.9794 

0.9489 

0.868 

3.5 

0.99975 

0.99965 

0.99965 

0.99941 

0.99905 

0.99799 

0.99581 

0. 9888 

0.9728 

4.0 

0.99993 

0.99995 

0.99995 

0. 99987 

0.  99977 

0.99965 

0. 99929 

0.  99853 

0. 99627 

4.5 

- 

- 

- 

0.99995 

0. 99995 

0.99995 

0.  99991 

0.  99983 

0.99967 
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Table  1(1).  The  Probability  Distribution  of  PA(yo>  z,  F  =  1),  Where  yQ  is  the 

(1/1 0)th  Maximum  of  the  Area  (A)  Given  by  n/A=  r.  2Z.  The  sample  size  was 
25.  000  maps 


Z 

0 

1 

2 

3 

4 

5 

6 

7 

8 

-3.5 

0.00015 

0.00015 

0.00017 

0.00001 

- 

- 

- 

- 

- 

-3.0 

0.00097 

0. 00115 

0.00103 

0.00069 

0.00053 

0.00009 

- 

- 

- 

-2.5 

0.00501 

0.00519 

0.  00475 

0.00329 

0.00221 

0.00029 

0.00003 

- 

- 

-2.0 

0.0222 

0.0209 

0.0185 

0.0142 

0.00929 

0.00383 

0.00039 

- 

- 

-1.5 

0.0628 

0.0612 

0.0549 

0.0460 

0.0328 

0.0166 

0.00331 

- 

- 

-1.0 

0.  154 

0.  149 

0.  138 

0.  123 

0.0924 

0.0569 

0.0190 

0.00067 

- 

-0.5 

0.302 

0.294 

0.277 

0.  251 

0.  206 

0. 144 

0.0697 

0.00787 

- 

0 

0.491 

0.482 

0.465 

0.430 

0.376 

0.297 

0.  193 

0.0616 

0.00035 

0.5 

0.680 

0.673 

0.655 

0.626 

0.578 

0.499 

0.389 

0.223 

0.0239 

1.0 

0.832 

0.826 

0.814 

0.792 

0.757 

0.699 

0.616 

0.499 

0.288 

1.5 

0.9319 

0.9267 

0.9204 

0. 9080 

0.889 

0.855 

0.814 

0.771 

0.768 

2.0 

0.9758 

0. 9737 

0. 9712 

0.9666 

0.9597 

0.  946  2 

0.9334 

0.9300 

0.9680 

2.5 

0.99431 

0.99341 

0.99241 

0. 99057 

0.  987  9 

0.9844 

0.9816 

0.9872 

0.  99827 

3.0 

0.99877 

0.99857 

0.99835 

0.99769 

0.99713 

0.99609 

0.  99603 

0.99817 

0.99997 

3.5 

0.99975 

0.99967 

0.99969 

0. 99965 

0.99929 

0.99941 

0.99945 

0.99985 

- 

4.0 

0.99993 

0.99995 

0. 99995 

0. 99995 

0.99991 

0.  99989 

0.  99997 

0.99995 

- 

4.5 

- 

- 

- 

- 

0. 99995 

- 

- 

- 

- 
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Table  1(2).  The  Probability  Distribution  of  PA  (yQ.  z,  F=  2).  where  yQ  is  the 
(2/l0)ths  Maximum  of  the  Area  (A)  Given  by  \lA=  2zr 


g=-.r  "  -  -  -  -■  -  r  : -  - 

2 

^o 

0 

1 

2 

3 

4 

5 

6 

7 

8 

-3.5 

0.00019 

0.00015 

0.00017 

0.00009 

0.00001 

- 

- 

- 

- 

-3.  0 

0.00107 

0.00115 

0.00111 

0.00087 

0.00077 

0.00005 

- 

- 

- 

-2.5 

0. 005 1 1 

0.00537 

0.  004  8  9 

0. 00385 

0.00275 

0.00091 

0.  00005 

- 

- 

-2.0 

0.0224 

0.0218 

0.0196 

0.0168 

0. 0121 

0.00643 

0.00103 

- 

- 

-1.5 

0.0634 

0.0630 

0.057  5 

0.0518 

0.0408 

0.0246 

0.00743 

0. 00035 

- 

-1.0 

0.  155 

0.  152 

0.  143 

0.  135 

0. 1101 

0.  0772 

0.0372 

0.00387 

- 

-0.5 

0.305 

0.299 

0.287 

0.267 

0.234 

0.  184 

0.  118 

0.00318 

0. 00007 

0 

0.493 

0.487 

0.476 

0.451 

0.413 

0.356 

0.279 

0.  150 

0.00861 

0.5 

0.683 

0.678 

0.665 

0.645 

0.612 

0.563 

0.500 

0.395 

0. 183 

1.0 

0.834 

0.830 

0.821 

0.806 

0.7  84 

0.755 

0.720 

0.689 

0.703 

1.5 

0.  9330 

0. 9285 

0.9243 

0.9157 

0.  9051 

0.888 

0.882 

0.891 

0.9656 

2.0 

0.9761 

0.9748 

0.9726 

0. 9704 

0.9664 

0.9617 

0.9620 

0.9774 

0.99865 

2.5 

0.  99445 

0.99383 

0.99287 

0.99189 

0.99027 

0.9896 

0.99203 

0.99683 

- 

3.0 

0.99879 

0.99863 

0.99845 

0.99821 

0.99769 

0.99743 

0.99825 

0.99973 

- 

3.5 

0.99979 

0.99975 

0.99973 

0. 99977 

0. 99955 

0.99965 

0. 99995 

0.99995 

- 

4.0 

0.99993 

0. 99995 

- 

- 

0.99993 

- 

- 

- 

- 

4.5 

- 

" 

- 

- 

“ 

- 

“ 
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Table  1(3).  The  Probability  Distribution  of  PA(yQ.  z.  F*  3).  Where  yQ  is  the 
(3/10)ths  Maximum  of  the  Area  (A)  Given  by  *^A  *  2Z.  r 


z 

^o 

0 

1 

2 

3 

4 

5 

6 

7 

8 

-3.  5 

0.  00019 

0.00015 

0.  00021 

0. 00009 

0. 00005 

- 

- 

- 

- 

-3.0 

0.00111 

0.  00121 

0.00113 

0. 00093 

0. 00089 

0.00023 

- 

- 

- 

-2.5 

0. 00527 

0.  00553 

0.00517 

0.00443 

0.00335 

0.00153 

0.00019 

- 

- 

-2.0 

0. 0227 

0.0223 

0,  0206 

0.0186 

0.0150 

0.00919 

0.00235 

0.  00003 

- 

-1.5 

0.  0641 

0. 0334 

0. 0600 

0.  0569 

0.  0490 

0. 0329 

0.0144 

0.00107 

- 

-1.0 

0.  158 

0.  156 

0.  149 

0.  143 

0.  126 

0. 0984 

0.  0589 

0.0114 

- 

-0.  5 

0.  307 

0.  304 

0.  295 

0.  281 

0.  258 

0.  221 

0.  165 

0.0734 

0.00 101 

0 

0.495 

0.  492 

0.  484 

0.  468 

0.  444 

0.407 

0.  355 

0.  258 

0.0600 

0.  5 

0.  685 

0.  682 

0.  673 

0.  662 

0.  642 

0.  615 

0.  585 

0.546 

0.  493 

1.  0 

0.  836 

0.  833 

0.  828 

0.  817 

0.  806 

0.794 

0.789 

0.  810 

0.  9225 

1.  5 

0. 9333 

0. 9307 

0.  9276 

0. 9230 

0.  9168 

0.  9122 

0.9200 

0. 9500 

0. 99665 

2.0 

0.  9763 

0.  9754 

0.  9739 

0. 9739 

0.9720 

0.  9723 

0. 9783 

0. 99285 

0. 99991 

2.5 

0. 99455 

0. 99409 

0.  99317 

0.  99281 

0. 99213 

0.99245 

0.99575 

0. 99907 

- 

3.0 

0. 99879 

0.  99865 

0. 99853 

0. 99839 

0.  99813 

0. 99837 

0. 99923 

0. 99995 

- 

3.  5 

0. 99979 

0.  99975 

0.  9997  5 

0. 99977 

0. 9997  1 

0.  99979 

0. 99995 

- 

- 

4.0 

0. 99993 

0. 99995 

- 

- 

0. 99995 

- 

- 

- 

- 

4.  5 

- 

“ 

** 

“ 

Table  1(4).  The  Probability  Distribution  of  PA  (yQ,  z,  F  =  4),  Where  y  is  the 
(4/  10)ths  Maximum  of  the  Area  (A) 


— 

Z 

m 

0 

1 

2 

3 

4 

5 

6 

7 

8 

-3.  5 

0.00019 

0.00017 

0.00021 

0.00009 

0.00013 

0. 00003 

- 

- 

- 

-3.0 

0.00111 

0.00123 

0.00117 

0.  00107 

0. 00097 

0.  00029 

- 

- 

- 

-2.5 

0.00533 

0. 00563 

0. 00553 

0.  00491 

0.00417 

0.00257 

0.  00049 

- 

- 

-2.0 

0. 0229 

0.0228 

0.0218 

0.0204 

0.0180 

0.0123 

0.00437 

0. 00025 

- 

-1.  5 

0. 0645 

0.  0655 

0.0632 

0.0618 

0. 0559 

0.0440 

0.0238 

0.00307 

- 

-1.0 

0.  159 

0.  159 

0.  155 

0.  152 

0.  142 

0.  120 

0.  0848 

0. 0278 

0.  00003 

-0.  5 

0.  309 

0.  308 

0.  303 

0.  295 

0.  281 

0.  260 

0.  222 

0.  136 

0.  0118 

0 

0.  498 

0.  496 

0.  492 

0.  483 

0.  472 

0.  455 

0.  429 

0.  377 

0.  225 

0.  5 

0.  687 

0.  686 

0.  682 

0.  675 

0.  671 

0.  660 

0.  656 

0.  673 

0.  781 

1.0 

0.  838 

0.  836 

0.  834 

0.  829 

0.  825 

0.  824 

0.  843 

0.  888 

0. 9853 

1.  5 

0. 9340 

0.  9319 

0.  9311 

0. 9285 

0. 9259 

0. 9308 

0. 9444 

0.  9791 

0.  99995 

2.0 

0.  9765 

0.  9761 

0.  9759 

0.  9760 

0.  9759 

0.  9788 

0.  9865 

0. 99769 

- 

2.  5 

0. 99455 

0.  99417 

0. 99353 

0. 99359 

0. 99393 

0. 99477 

0. 99787 

0.  99975 

- 

3.0 

0. 99879 

0.  99867 

0. 99857 

0. 99849 

0.  99859 

0.  99891 

0. 99975 

- 

- 

3.  5 

0. 99979 

0. 99977 

0.9997  5 

0. 99977 

0.99971 

0.  99989 

- 

- 

- 

4.0 

0. 99993 

0. 99995 

- 

- 

0. 99995 

- 

- 

- 

- 

4.  5 

- 

_ 

_ 

_ 

_ 

_ 

_ 

_ 
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Table  1(5).  The  Probability  Distribution  of  PA(yQ,  z,  F=5),  Where  yQ  is  the 
(5/10)ths  Maximum  of  the  Area  (A)  Corresponding  to  Z  =  0(1)8,  A=  256 


Z 

y0 

0 

1 

2 

3 

4 

5 

6 

7 

8 

-3.5 

0.00019 

0.00017 

0.00023 

0.00011 

0. 00021 

0. 00009 

- 

- 

- 

-3.0 

0. 00119 

0.00131 

0.  00121 

0. 00141 

0. 00113 

0. 00053 

0. 00007 

- 

- 

-2.5 

0.00539 

0.00571 

0.006  1  3 

0.00569 

0.004  8  7 

0.00421 

0.00115 

0.00001 

- 

-2.0 

0.0231 

0.0234 

0.0227 

0. 0225 

0.0208 

0.0162 

0.00801 

0.00085 

- 

-1.5 

0.0654 

0.0671 

0. 0656 

0.0666 

0.0642 

0.0548 

0.0378 

0.00845 

- 

-1.0 

0.  160 

0.  162 

0.  160 

0.  161 

0.  157 

0.  146 

0.  117 

0.0591 

0. 00147 

-0.5 

0.311 

0.311 

0.310 

0.309 

0.305 

0.  298 

0.  280 

0.  222 

0. 0606 

0 

0.500 

0.500 

0.  500 

0.500 

0.500 

0.500 

0.  500 

0.500 

0.500 

0.5 

0.689 

0.689 

0.690 

0.691 

0.695 

0.702 

0.  7  20 

0.778 

0. 9394 

1.0 

0.840 

0.838 

0.  840 

0.839 

0.843 

0.854 

0.  883 

0. 9409 

0. 99851 

1.5 

0.9346 

0.9329 

0.9343 

0.9333 

0.9358 

0. 9452 

0.9622 

0.  99153 

- 

2.0 

0.9769 

0.9766 

0.9773 

0.9774 

0. 9791 

0. 9838 

0.  99197 

0.  99913 

- 

2.5 

0.99459 

0.99427 

0.99385 

0.99429 

0.99511 

0.99577 

0.99883 

0.99997 

- 

3.0 

0.99879 

0. 99867 

0.99877 

0.99857 

0.99885 

0. 99945 

0.99991 

- 

- 

3.5 

0.99979 

0.99981 

0.99975 

0.99987 

0.  99977 

0.  99991 

- 

- 

- 

4.0 

0. 99993 

- 

- 

- 

- 

- 

- 

- 

- 

4.5 

• 

- 

- 
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Figure  4(0).  For  the  Graphical  Solution  of  P  (yQ,  z,  F),  the  Probability  That  an 
END  (y)  Will  Exceed  (y0)  Only  in  (F/10)  or  Less  of  the  Area  (A)  When  Scale 
Distance,  r=  1.  (This  chart  is  for  F  =  0. )  The  horizontal  scale  is  for  the 
standardized  measure  (z)  of  the  area  [Eq.  (21)] 


f.  y  \  the  Probability  ^bat  an 

Figure  4(4).  For  the  Graphical  Solution  of  *>'Area  (A>  When  Scale 

ENgD  (y)  Will  Exceed  (y  \£rtie'tofl  =  4. )  The  horizontal  scale  is  for  the 
Standardized  Las^re  U)  of  the  area  (Eq.  (21)] 
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Figure  4(5).  For  the  Graphical  Solution  of  P  (yQ,  z,  F)f  the  Probability  That  an 
END  (y)  Will  Exceed  (yQ)  Only  in  (F/10)  or  Less  of  the  Area  (A)  When  Scale 
Distance,  r=  1,  (This  chart  is  for  F  =  5, )  The  horizontal  scale  is  for  the 
standardized  measure  (z)  of  the  area  [Eq.  (21)] 
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2.4.3  CLOUD-COVER  DISTRIBUTION  BY 
GRAPHICAL  SOLUTION 

Begin  with:  mean  cloud  cover: 

scale  distance:  r  km 

2 

area  of  floor  space:  A  km  . 

(For  the  sky  dome,  as  the  surface  observer  sees  it,  the  floor  space  is  assumed  to 

2 

have  an  area  of  A  =  2424  km  equivalent  to  a  circular  area  with  a  15-nm  radius. ) 
Find 

z  =  ln(sl  A/r)/fn  2  . 

Find  the  cumulative  probability 

P  =  1  -  P '  (25) 

o  o 

which  is  also  the  single-point  probability  of  no  cloud.  Find  its  END,  y  .  On  each 

chart,  for  each  F  =  0(1)5,  the  intersection  of  the  y  -curve  with  the  ordinate  line 

o 

through  z,  interpolated  as  necessary,  gives  the  reading  of  the  probability, 

PA(yQ.  z,  F). 

For  F  =  6(1)10  we  take  advantage  of  mirror-image  symmetry.  Enter 
Figure  4(4)  down  to  Figure  4(0)  respectively,  and  find  the  intersection  of  the  curve 
of  (-yQ)  with  the  ordinate  line  through  z.  Read  the  probability  PA  (yQ,  z,  F)  on 
the  right-hand  scale. 

Or,  for  F  =  6(1)10  we  can  approach  the  problem  as  follows:  look  for  the  solu¬ 
tion  for  PA  (-yQ,  A,  (10-F),  r)  on  the  (lO-F)-chart  at  the  intersection  of  the 
( -yQ)-curve  with  the  z-ordinate.  Then, 

PA  (yQ,  A,  F,  r)  =  1  -  PA  (-yQ.  A,  (10-F),  r)  .  (26) 

Example  1:  Suppose  we  wish  to  find  the  probability  that  the  single-point  lower 

10  percent  of  visibility  is  to  be  exceeded  in  no  more  than  (2/10)ths  of  a  region  of 
2 

area  250  km  .  The  scale  parameter  is  given  (say,  r  =  2  km).  Then 
z  =  In  (VA/r)/fn  2  =  2.  98  , 

PQ  =  0.1,  for  which  yQ  -1.  28, 

and 

F  =  2  . 
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Entering  Figure  4(2),  (the  dashed  curves  and  lines  are  for  the  next  example)  and 
following  the  curve  for  yQ  =  -1.  28  to  the  intersection  with  z  =  2.  98,  we  obtain,  on 
the  left-hand  scale,  the  probability 

PA(yo,  A.  F.  r)  =  0.  082  . 

That  is,  the  answer  is  8.  2  percent,  only  slightly  lower  than  the  basic  single-point 

o 

probability.  If  this  visibility  were  not  to  be  exceeded  anywhere  in  the  250  km 
area  the  answer  [from  Figure  4(0)]  would  have  been  6  percent. 

1 3 

Example  2:  For  Hanscom  AFB,  MA  the  RUSSWO  gives  a  mean  sky  cover  of 
0.66  in  January,  noontime,  assume  for  this  example,  that  the  scale  distance 
r  =  1 .  9  km  then, 

P'  =  0.66,  P  =  0.  34,  y  =  -0.41  . 
o  o  J  o 

2 

For  the  sky  cover,  the  area  A  is  taken  to  be  2424  km  which,  together  with 
r  =  1.  9  km,  gives 

z  =  in  (V2424/1.9)/fn2  =  4.70  .  (27) 

By  finding  the  curve  y^  =  -.041  in  each  Figure  4(0)  to  Figure  4(5)  successively, 
and  following  the  curve  until  z  =  4.  70,  the  probabilities  of  all  clear,  1  / 1 0  cover, 

- ,  (5/10)ths  cover  were  found,  as  shown  in  the  column  for  estimates  by  graph 

(Table  2).  By  entering  Figures  4(4)  down  to  Figure  4(0)  successively  at 
yo  =  0.41  and  following  each  curve  until  z  =  4.70,  the  probabilities  of  6/10  to  full 
overcast  were  found  on  the  righthand  scale,  as  shown  in  the  rest  of  Table  2.  The 
RUSSWO  figures  are  shown  for  comparison.  The  broken  lines  for  this  example 
have  been  plotted  in  Figures  4(0)  to  Figure  4(5). 


RUSSWO  -  Revised  Uniform  Summaries  of  Surface  Weather  Observations. 
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Table  2.  Frequency  Distribution  of  Sky  Cover  at  Hanscom  AFB,  MA, 
January,  Noontime,  as  Given  in  the  RUSSWO,  and  as  Estimated  by  the 
3D-BSW  Model.  The  mean  sky  cover  is  P  =  0.66,  scale  distance, 

o  2 

r  =  1.9  km.  The  sky-dome  area  is  assumed  to  be  A  =  2424  km  , 
giving  z  =  4.70 


! - 

Sky 

Cover 

F 

RUSSWO 

3D- 

BSW  Estimates 

By  Graph 

By  Computer  (Appendix  B) 

Freq. 

Freq. 

Freq. 

Clear 

0 

0.  155 

4(0) 

0.  140 

0.  127 

1/10 

1 

0.  198 

4(1) 

0.  190 

0.  191 

§  2/10 

2 

0.  247 

4(2) 

0.  230 

0.  230 

5  3/10 

3 

0.  288 

4(3) 

0.  262 

0.  266 

5  4/10 

4 

0.  329 

4(4) 

0.  297 

0.  296 

5  5/10 

5 

0.  360 

4(5) 

0.  325 

0.  330 

S  6/10 

6 

0.  387 

4(4) 

0.  375 

0.  366 

5  7/10 

7 

0.431 

4(3) 

0.400 

0.401 

5  8/10 

8 

0.491 

4(2) 

0.440 

0.446 

5  9/10 

9 

0.537 

4(1) 

0.51 

0.499 

5  10/10 

10 

4(0) 

0.61 

0.605 

Overcast 

0.462 

0.  39 

0.  395 

2.4.4  CLOUD  COVER  DISTRIBUTION  BY  ALGORITHM 
AND  COMPUTER  SOLUTION 

4 

Algorithms,  previously  published  with  the  2D-BSW  model,  were  found  to 
approximate  the  answers  of  the  3D-BSW  model  with  rmse  =  0.005,  bias  -0.003, 

If  limited  to  z  5  7,  the  bias  is  as  low  as  -0,  001.  This  high  degree  of  accuracy 
permits  us  to  use  the  previous  algorithms  (Appendix  B)  for  the  probability  of  areal 
coverage,  PA  (y  ,  A,  F,  r). 

Example  2  (Contd, );  Table  2  shows,  in  the  last  column,  the  results  of  solution  by 
computer  (following  the  technique  described  in  Appendix  B).  The  values  are  close 
to  the  values  bv  graphical  solution.  (The  graphical  values  might  be  read  slightly 
differently  by  another  user  of  the  graphs.  ) 
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The  difference  in  the  estimates  (by  graphs  or  by  algorithms)  from  the 
RUSSWO  frequencies  can  be  party  explained  by  the  fact  that  the  RUSSWO  categories 
include  up  to  5  percent  additional  cover,  except  for  the  overcast  that  includes 
breaks  of  5  percent.  The  model,  on  the  other  hand,  gives  probability  of 
strictly  clear,  and  strictly  overcast  without  breaks.  For  the  same  reason  each 
partial  cover  should  appear  slightly  less,  in  cumulative  frequency,  than  in 
the  RUSSWO,  This  effect  is  clearly  illustrated  in  Figure  5.  The  latter  shows 


CLOUD  COVER  (F)-TENTHS 


Figure  5.  A  Plot  of  the  Cumulative  Probability,  PR(S  F),  vs  Cloud 
Cover  (0  ^  F  £  10)  in  Tenths.  The  X's  are  plotted  values  from  the 
RUSSWO  for  Hanscom  AFB,  January,  noontime,  with  mean  0.66. 
With  scale  distance  determined  at  r  =  1.  9  km,  the  3D-BSY,  model 
provided  the  three  curves  for  areas  (10  km)^,  2424  km^  (sky  cover) 
and  (100  km)2 


\ 
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the  plot  of  the  RUSSWO  sky-cover  frequencies  (at  the  X's),  the  3D-BSW  model  prob- 

2 

ability  curve  for  the  sky  cover  of  area  2424  km  ,  and,  in  addition,  the  model  curves 

2  4  2 

for  cloud  cover  over  areas  of  100  km  and  10  km  .  The  additional  two  curves  are 
estimates  without  benefit  of  verification. 

25  Line  Coverage  (PL) 

The  probability  distributions  of  conditions  along  a  line,  or  fraction  of  a  line, 
were  obtained,  together  with  algorithms,  for  solution  by  the  previous  2D-BSW  model. 
This  exercise  is  not  repeated  for  the  3D-BSW  model,  since  a  more  important  goal 
takes  precedence,  as  described  in  the  following  section.  For  an  algorithm  solution 
of  PL(yn,  s',  F,  r),  (the  probability  that  yQ  will  be  exceeded  in  no  more  than 
(1  /l0)ths  of  the  lineof  length  (s'))  based  on  results  with  the  2D-BSW  model  see 
Appendix  C. 

2.6  Line  Intervals  (PI) 

The  goal  of  this  section  is  to  estimate  the  probability,  PI(X,  s',  T,  r),  that  a 
condition  (X)  is  the  minimized  maximum  in  a  linear  interval,  s' (km),  somewhere 
along  a  longer  line  of  travel  of  length  T(km).  The  given  parameter  is  either  the 
scale  distance  r(km)  or  the  wavelength  A(km)  where  A  =  256  r,  in  the  3D-BSW  model. 
In  terms  of  s',  T,  and  r,  the  following  standardized  terms  are  obtained: 

s  =  s'/r 

z  =  fn  s/fn  2  (28) 

(j)  —  C n  (T/s')/fn  2  . 

With  r  =  1  or  s'  =  s,  the  Monte  Carlo  exercise,  previously  described,  was  used  to 
produce  25,  000  maps  (250,  000  would  have  been  better  but  prohibitive),  each  with 
24  wave  formations  that  produced  a  field  of  END-  or  y-values  throughout  the  3D 
space.  For  this  exercise  the  y-values  were  determined  along  a  line,  arbitrarily 
chosen  to  be  the  V-axis  (Figure  2).  The  line  interval  (s)  was  made  equal  to  2Z  for 
z  -1(1)8.  The  overall  line  of  travel  (T)  was  made  equal  to  2Z+U>,  u>=  0(l)(10-z),  so 
that  the  length  of  travel  (T)  varied  from  a  minimum  equal  to  one-half  of  the  single¬ 
unit  interval  to  a  maximum  of  1024  units  of  distance.  The  y-values  were  generated 
along  the  line,  equally  spaced  at  small  intervals: 

6s  =  s/N  .  (29) 


The  number  N  was  arbitrarily  set  at  20,  so  that  y-values  were  determined  at  21 
points  along  the  length  of  (s). 
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The  interval  (s)  was  moved  along  the  line  (T)  to  find  the  lowest  maximum  in 
that  interval.  This  was  done  25,000  times,  yielding  a  frequency  distribution  of 
y(minmax)  for  each  combination  of  s  and  T,  or  z  and  ux 

Tables  3(-l)  to  3(8)  show  the  results  for  the  frequencies  that  approximate  the 
probability  distribution,  PI  (y  ,  z,  o»),  applicable  to  s'  and  T  when  the  scale  distance 
(r)  is  given,  thus: 

s'  =  r.  2Z  for  z  =  -1(1)8 

(30) 

T  =  r.  2Z  +  U*  for  u>  =  0(l)(10-z) 

2.  6.  1  BY  GRAPHICAL  SOLUTION 

Figures  6(  —  1 )  to  6(8)  show  the  same  results,  except  that  the  curves  have  been 
smoothed  as  deemed  desirable,  especially  for  the  very  high  and  very  low  probabilities. 
Each  figure  is  for  a  single  value  of  z.  In  practice  there  would  need  to  be  a  system 
of  interpolation  between  the  values  in  the  tables  or  figures. 

2.  6.  2  BY  ALGORITHM  AND  COMPUTER  SOLUTION 

While  the  figures  could  be  used  directly,  it  is  much  more  desirable,  in  fact 

imperative,  to  compose  algorithms  for  the  estimates  of  probabilities.  This  was  done, 

with  results  as  contained  in  Appendix  D.  The  basic  equation  is  of  the  form: 

2 

y  =  A  (z,  m)  +  B(z,  cu)  .  yQ+  C  (z,  cu)  .  y  .  After  the  y  -  curves  were  drawn  of  P  vs  uo, 
for  each  Z  =  -1(1)8  they  were  extrapolated  beyond  (10-Z)  to  (12-Z),  by  eye,  to  give 
a  useful  extension  of  the  model's  application.  The  additional  values  of  A(i,  j), 

B(i,  j),  C  (i,  j)  are  included  in  the  tables  (Appendix  D).  Example:  Suppose  that  a 
cloud-free  interval  of  1  km  is  sought  in  a  flight  track  of  length  of  50  km,  when  the 

7 

sky  cover  (f)  is  observed  to  be  (9/10)ths.  The  Stanford  Research  Institute  (SRI) 
formula  for  a  cloud-free  line  of  sight  (CFLOS)  overhead  gives 

P  =  1  -  f  (1+  3f) / 4  =  0.  168  . 
n 

Accepting  this  as  the  single-point  probability  (Pq),  then  yQ  =  -0.  96. 

Suppose,  further,  that  the  parameter  value  is  r  -  0.  109  km,  or  A  =  27.6  km  (see 
below).  Using  Appendix  D,  with  this  value  of  r,  together  with  s'  =  1  km,  T  =  50  km, 
for  which  we  obtain  z  =  3.21,  u=  5.64,  the  solution  to  the  algorithms  gives 

PI  (s'  =  1  km,  T  =  50  km)  =  0.  84  . 

That  is,  there  is  an  84  percent  probability  of  having  a  1-km  window  of  observation 
through  the  clouds  on  a  50-km  track,  when  the  mean  cloud  cover  is  (9/10)ths. 

7.  Allen,  J.  H. ,  and  Malick,  J.  D.  (1983)  The  frequency  of  cloud-free  viewing 
intervals.  Twenty-first  Aerospace  Science  Meeting,  10-13  .January  1983, 

Reno,  NV. ,  Copyright  AIAA  Inc. 
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Table  3 < - 1 ) .  Estimate  of  the  Probability  Distribution,  Pr(y  ,  z  =  -1,  u>),  of  the  Minimized  Maximum  (y  )  i 


0.999767  3.5  0.  99977  0.999' 


Table  3(0).  Estimates  of  the  Probability  Distribution,  PI  (y  ,  0,  u>),  of  the  Minimized  Maximum  (y 
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Table  3(2).  Estimates  of  the  Probability  Distribution.  PI  (yQ,  z  =  2,  v),  of  the 
Minimized  Maximum  (y  )  in  a  Linear  Interval  (S  =  2  =4)  Over  a  Total  Line  of 
Travel  (T  =  2Z+U)  =4  to  1024  units) 


lea 

j) 

0 

1 

9 

3 

4 

5 

6 

7 

8 

-4.0 

- 

- 

- 

- 

0.00014 

0.00018 

0. 00050 

0.00090 

0.00170 

-3.5 

- 

0.00018 

0. 00026 

0. 00054 

0. 00098 

0. 00150 

0. 00254 

0. 00514 

0. 00986 

o 

• 

CO 

0.00114 

0.00142 

0.  00202 

0.00278 

0.  00494 

0.00846 

0.0166 

0.0305 

0.0573 

-2.5 

0.  00530 

0.  007  38 

0.00926 

0.  0132 

0. 0197 

0.  0345 

0.0641 

0.  1 19 

0.217 

-2.0 

0.0211 

0.0238 

0.0300 

0.0424 

0. 0649 

0.  107 

0.  193 

0.  337 

0.545 

-1.5 

0.0609 

0.0672 

0.0855 

0.  1 15 

0.  170 

0.  259 

0.4  24 

0.642 

0.857 

-1.0 

0.  149 

0.  167 

0.  199 

0.  252 

0.336 

0.478 

0.  6  95 

0.886 

0. 9835 

-0.5 

0.292 

0.  318 

0.  361 

0.433 

0.544 

0.705 

0.892 

0. 9832 

0.  9950 

0 

0.479 

0.504 

0.554 

0.6  34 

0.750 

0.878 

0. 9753 

0. 99846 

0.  99998 

0.5 

0.672 

0.6  95 

0.742 

0.807 

0.885 

0. 9615 

0.996 94 

0. 99994 

- 

1.0 

0.826 

0.844 

0.877 

0.9153 

0. 9601 

0.  99302 

0. 99998 

0. 99998 

- 

1.5 

0.9243 

0.  9339 

0. 9504 

0. 9706 

0.  9892 

0.  998  98 

- 

- 

- 

2.0 

0.9744 

0.9784 

0. 9835 

0.99246 

0.99758 

0. 99982 

- 

- 

! 

2.5 

0.99286 

0.  99474 

0.  99650 

0.  99890 

0.99970 

0. 99998 

- 

- 

1 

! 

3.0 

0. 99838 

0.  99890 

0.  99918 

0. 99970 

0.99998 

- 

- 

- 

- 

3.5 

0.99970 

0.99970 

0.99974 

0. 99994 

- 

- 

- 

- 

i 

4.0 

0.99998 

0. 99998 

0. 99998 

0. 99998 

- 

- 

- 

- 

- 
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Table  3(3).  Estimates  of  the  Probability  Distribution.  PI  (y  ,  z  =  3,  u>),  of  the 
Minimized  Maximum  (y^)  in  a  Linear  Interval  (S  =  2  =8)  Over  a  Total  Line  of 
T ravel  (T  =  2Z  +  W=  8  to  1024  units) 


0.00006  0.00006  0.00014  0.00014  0.00050 

0.00078  0.00086  0.00154  0.00354  0.00646 

0.00462  0.007  10  0.0125  0.0221  0.0423 

0.0152  0.0267  0.0508  0.0963  0.178 

0.0533  0.0902  0.162  0.289  0.484 

0.  141  0.224  0.374  0.585  0.823 

0.291  0.  433  0.645  0.855  0.9772 

0.494  0.665  0.861  0.9740  0.99906 

0.  703  0.  848  0.  9672  0.  99782  0.  99998 


-3.  5 

0.00002 

0.00002 

0. 00054 

0. 00078 

0.  00086 

0.  00154 

0.  00354 

-3.  0 

0.00134 

0.  00190 

0.  00306 

0. 00462 

0.  007  10 

0.0125 

0.0221 

-2.  5 

0. 00414 

0. 00622 

0. 00922 

0.  0152 

0.  0267 

0. 0508 

0. 0963 

-2.  0 

0.  0  184 

0.  0240 

0.  0340 

0. 0533 

0.0902 

0.  162 

0.  289 

-1.  5 

0. 0535 

0.  067  1 

0. 0945 

0.  141 

0.  224 

0.  374 

0.  585 

-1.  0 

0.  131 

0.  165 

0.  210 

0.  291 

0.  433 

0.  645 

0.  855 

-0.  5 

0.  264 

0.  309 

0.  383 

0.  494 

0.  665 

0.  861 

0. 9740 

0 

0.  451 

0.  498 

0.  582 

0.703 

0.  848 

0.  9672 

0. 99782 

0.  5 

0.  643 

0.  689 

0.  765 

0.  857 

0. 9463 

0. 99530 

0.  99982 

1.0 

0.  802 

0.  836 

0.  887 

0. 9443 

0. 9889 

0. 99966 

0. 99998 

1.  5 

0. 9132 

0.  9313 

0.  959  1 

0. 9833 

0.  99778 

0.  99998 

- 

2.  0 

0. 9707 

0. 9774 

0. 9895 

0. 99734 

0. 99990 

- 

- 

2.  5 

0.  99150 

0.  9937  8 

0.  99790 

0.  99970 

0.  99998 

- 

- 

3.  0 

0.  99798 

0.  99870 

0.  99970 

0.  99982 

- 

- 

- 

3.  5 

0.  99982 

0. 99982 

0.  99998 

0.  99998 

- 

- 

- 

4.  0 

0.  99998 

0. 99998 

_ 

_ 

_ 

Table  3(4).  Estimates  of  the  Probability  Distribution,  PI  (y  ,  z  =  4,  uj),  of  the 
Minimized  Maximum  (y  )  in  a  Linear  Interval  (S=  2Z=  16)  Over  a  Total  Line  of 
Travel  (T  =  2Z  16  to  1024  units) 


It) 

0 

1 

2 

3 

4 

5 

6 

-4.0 

- 

- 

0.  00006 

0. 00006 

0.  00014 

0.  00014 

0. 00042 

-3.  5 

- 

0.  000  10 

0.  00018 

0.  00026 

0. 00066 

0.  00154 

0. 00342 

-3.  0 

0. 00094 

0.  00186 

0. 00250 

0. 00454 

0. 00746 

0.  0128 

0. 0247 

-2.5 

0.00298 

0. 00530 

0.  00878 

0.  0  168 

0.  0327 

0.  0611 

0.  114 

-2.  0 

0.  0135 

0.0213 

0. 0355 

0.  0624 

0.  118 

0.  212 

0.  377 

-1.  5 

0. 0409 

0.  0605 

0. 0977 

0.  162 

0.  289 

0.  485 

0.  731 

-1.0 

0.  105 

0.  149 

0.  224 

0.  351 

0.  553 

0.  782 

0. 9507 

-0.  5 

0.  222 

0.292 

0.  403 

0.  57  6 

0.798 

0. 9502 

0.  99718 

0 

0.  402 

0.  481 

0.  611 

0.  783 

0.  9416 

0. 99474 

0. 99998 

0.  5 

0.  601 

0.  678 

0.  797 

0.  9165 

0.  9899 

0. 99966 

- 

1.0 

0.  773 

0.  831 

0.  9133 

0. 9790 

0.  99906 

0.  99998 

- 

1.  5 

0.  897 

0. 9286 

0.  97  12 

0.  99586 

0.  99998 

- 

- 

2.0 

0. 9633 

0.  9780 

0.  99434 

0.  99982 

- 

- 

- 

2.  5 

0. 9897 

0. 99422 

0.  99914 

0.  99998 

- 

- 

- 

3.0 

0.  997  54 

0.  99902 

0.  9997  8 

- 

- 

- 

- 

3.  5 

0. 99978 

0.  99994 

0.  99998 

- 

- 

- 

- 

4.0 

0.  99982 

0.  99998 

- 

_ 

_ 

__ 

39 


40 


Tables  3(7)  and  3(8).  Estimates  of  the  Probability  Distribution, 
PI(yo;  z  =  7  or  8,  o>),  of  the  Minimized  Maximum  (y  )  in  a 

Linear  Interval  (S=2  =  128  or  256)  Over  a  Total  Line  of  Travel 

(T  =  2  <1,=  128  or  256  to  1024  units) 


■ - 

B 

■Mi 

cu(z 

=  7) 

w(z  =  8) 

IB 

0 

1 

2 

3 

0 

1 

2 

-4.0 

0.00006 

- 

- 

- 

_ 

_ 

-3.  5 

0. 00006 

- 

- 

- 

. 

-3.  0 

0.  00006 

- 

- 

_ 

1 

-  1 

-2.  5 

0. 00006 

- 

0.00002 

0.  00018 

- 

_ 

1 

-2.0 

0. 00038 

0.  00102 

0.  00174 

0. 00306 

- 

- 

0. 00006 

-1.  5 

0. 00098 

0. 00482 

0.0123 

0.  0260 

- 

- 

0. 00010 

-1.  0 

0. 00694 

0.  0256 

0.  0691 

0.  145 

0.  00046 

0.  00186 

0.  00366 

-0.  5 

0.  0327 

0.  103 

0.  248 

0.  452 

0.  00374 

0.  0146 

0. 0399 

0 

0.  116 

0.  27  6 

0.  552 

0.  824 

0.  0249 

0.  0901 

0.  211 

0.  5 

0.  'll  1 

0.  523 

0.  838 

0.  9810 

0.  109 

0.  297 

0.  581 

1.  0 

0.492 

0.  755 

0.  9687 

0.  99958 

0.  305 

0.  610 

0.  896 

1.  5 

0.720 

0.  9050 

0. 99686 

0.  99998 

0.  580 

0.  856 

0. 99074 

2.  0 

0.  880 

0.  9707 

0.  99978 

- 

0.  815 

0. 9648 

0.  999 ”7 

2.  5 

0. 9625 

0.  99318 

0. 99998 

- 

0.  9410 

0. 99222 

0.  99998 

3.  0 

0.  99162 

0.  99886 

- 

- 

0. 9858 

0. 99830 

_ 

3.  5 

0. 99934 

0. 99982 

- 

- 

0. 99822 

0.  99978 

_ 

4.0 

0. 99982 

0. 99994 

- 

0. 99950 

0. 99998 

- 
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Figure  6(-l).  The  Graphical  Solution  of  P  (yQ,  z,  w),  the  Probability 

That  (y  )  is  the  Minimized  Maximum  of  y  in  a  Line  Interval  s  (=  2Z) 

Throughout  a  Total  Length  of  Travel  (T  =  2Z+W),  Which  Become  Real 
Measures  When  Multiplied  by  the  Scale  Distance,  r(km). 

This  chart  is  for  z  =  -1:  T  ranges  from  1/2  to  1024  standardized  units 


Figure  6(0).  The  Graphical  Solution  of  P  (yQ,  z,  u>),  the  Probability 

That  (y  )  is  the  Minimized  Maximum  of  y  in  a  Line  Interval  s  (=  2Z) 

Throughout  a  Total  Length  of  Travel  (T=  2Z  +  W),  Which  Become  Real 
Measures  When  Multiplied  by  the  Scale  Distance,  r  (km). 

This  chart  is  for  z  =  0:  T  ranges  from  1  to  1024  units 


.000 1 


Kigure  6(1).  The  Graphical  Solution  of  P(yQ,  z.  u>),  the  Probability 

That  (y  )  is  the  Minimized  Maximum  of  y  in  a  Line  Interval  s  (=  2Z) 

Throughout  a  Total  Length  of  Travel  (T=  2Z+U>),  Which  Become  Real 
Measures  When  Multiplied  by  the  Scale  Distance,  r(km). 

This  chart  is  for  z=  1:  T  ranges  from  2  to  1024  units 
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.0001 


Figure  6(3).  The  Graphical  Solution  of  P  (yQ,  z,  «),  the  Probability 

That  (y  )  is  the  Minimized  Maximum  of  y  in  a  Line  Interval  s  (=  2Z) 

Throughout  a  Total  Length  of  Travel  (T  =  2Z+U>),  Which  Become  Real 
Measures  When  Multiplied  by  the  Scale  Distance,  r(km). 

This  chart  is  for  z=  3:  T  ranges  from  8  to  1024  units 
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Figure  6(4).  The  Graphical  Solution  of  P(yQ,  z»  w),  the  Probability 

That  (y  )  is  the  Minimized  Maximum  of  y  in  a  Line  Interval  s  (=  2Z) 

Throughout  a  Total  Length  of  Travel  (T  =  2Z+U)),  Which  Become  Real 
Measures  When  Multiplied  by  the  Scale  Distance,  r(km). 

This  chart  is  for  z  =  4 


Figure  6(5).  The  Graphical  Solution 
of  P(y  ,  z,  w),  the  Probability  That 
(y  )  is°the  Minimized  Maximum  of 

°  z 

y  in  a  Line  Interval  s  (=  2  )  Through¬ 
out  a  Total  Length  of  Travel 

(T  =  2Z+W),  Which  Become  Real 
Measures  When  Multiplied  by  the 
Scale  Distance,  r(km). 

This  chart  is  for  z  =  5 


Figure  6(6).  The  Graphical  Solution 

of  P(y  ,  z,  u>),  the  Probability  That 

(y  )  is°the  Minimized  Maximum  of 

°  z 

y  in  a  Line  Interval  s  (=  2  )  Through 

out  a  Total  Length  of  Travel 

(T=  2Z+UJ),  Which  Become  Real 
Measures  When  Multiplied  by  the 
Scale  Distance,  r(km). 

This  chart  is  for  z  =  6 


Figure  6(7).  The  Graphical  Solution 
°f  P(yQ.  z»  to),  the  Probability  That 
(yo)  is  the  Minimized  Maximum  of 

y  in  a  Line  Interval  s  (=  2Z)  Through¬ 
out  a  Total  Length  of  Travel 

(T=  2  Ui),  Which  Become  Real 
Measures  When  Multiplied  by  the 
Scale  Distance,  r(km). 

This  chart  is  for  z  =  7 


Figure  6(8).  The  Graphical  Solution 
P  (yD.  to),  the  Probability  That 
(yo)  is  the  Minimized  Maximum  of 

y  in  a  Line  Interval  s  (=  2Z)  Through 
out  a  Total  Length  of  Travel 

(T=  2Z+(V ),  Which  Become  Real 
Measures  When  Multiplied  bv  the 
Scale  Distance,  r(km). 

This  chart  is  for  z  =  8 


3.  THE  SRI  ALTERNATIVE  FOR  LINE  INTERVAL 


*7 

Allen  and  Malick  of  the  Stanford  Research  Institute  (SRI)  proposed  a  model  in 
which  the  mean  length  of  a  cloud-free  ihterval  (fi  km)  and  the  mean  length  of  a  cloudy 
interval  (A  km)  both  depend  on  the  probability  of  the  cloud-free  line-of-sight 
(PCFLOS).  If  P  denotes  PCFLOS  and  Q  =  1  -  P,  then  the  SRI  formula  provides 

H  (P)  =  (14P  +  1.  7) 

for  P  S  1/2  (31) 

A  (P)  =  ^  (14P  +  1.7) 

and 

M<  P>  =  -§04Q+  1.7) 

for  P  >  1/2 

A  (P)  =  (14Q+  1.  7)  . 

For  a  given  length  (s')  within  the  track  length  (T),  the  SRI  formula  becomes 

P's  1  -  jl -exp(-s-/M)jT/(M+X)  for  T  >  (M  +  A)  (32) 

or 

P^  =  P.exp(-s'/M)+Qjl-exp(-(T-s')/A)j.  exp(-s'  //i)  for  T  <  (/u  +  A)  . 

g 

Malick,  Allen  and  Zakanycz,  with  satellite  data  on  cloud-free  and  cloudy  intervals, 
in  the  month  of  Jime,  1978,  were  able  to  plot  a  frequency  of  clear  runs,  for  lengths 
between  1  and  11  km,  cloud  cover  such  that  single  point  cloud-free  probability  is 
0.27,  We  read  the  frequencies  of  run  lengths,  as  shown  (Table  4),  and  found  the 
best  estimates  of  the  scale  distance,  for  each  run  length  in  the  3D-BSW  model,  by 
using  the  model  in  successive  trial-and-error  steps.  With  the  resulting  geometric 
mean  value  of  the  scale  distance  (r-  0.  109  km),  the  3D-BSW  model  probabilities 
were  estimated,  to  compare  them  with  the  plotted  frequencies,  as  shown  (Table  4). 


Malick,  J.  D. ,  Allen,  J.  H. ,  and  Zakanycz,  S.  (1979)  Calibrated  analytical 
modeling  of  cloud-free  intervals.  Proceedings  of  Conference  SPIE. 

Vol.  195,  Atmospheric  Effects  on  Radiative  Transfer  (1979). 
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Table  4.  Readings,  From  an  SRI  Diagram,  of  the  Frequency  of 
Cloud-free  Run  Length,  When  the  Single -point  Probability  of  No 
Cloud  is  0.  27.  Included  also  are  estimates  of  corresponding 
scale  distance  in  the  3D-BSW  model,  and  the  latter's  probability 
estimates 


Run  Length 
(km) 

Frequency 

(SRI) 

Scale 
Distance 
r  (km  ) 

3D-BSW  ; 

Probability  1 

Estimate 
r  =  0.  109  km 

1 

0.  27 

0.  50 

0.  22 

2 

0.  25 

0.  60 

0.  18 

3 

0.  14 

0.  10 

0.  15 

4 

0.  12 

0.  11 

0.  12 

5 

0.  06 

0.07 

0.  10 

6 

0.  06 

0.08 

0.  086 

7 

0.  025 

0.  06 

0.075 

8 

0.  025 

0.  06 

0.  051 

9 

0.  025 

0.  07 

0.051 

10 

0.  02 

0.07 

0.  043 

11 

0.015 

0.  07 

0.  036 

Geometric  mean  =  0.  109 

The  greatest  departure  of  estimate  is  for  1-km  and  2-km  run  lengths.  Actually, 
since  the  single-point  probability  is  taken  to  be  0.  27,  the  1-km  cloud-free  run  should 
be  significantly  less,  but  was  not  recorded  that  way,  since  the  pixels  in  the  satellite 
pictures,  in  1-km  squares,  were  treated  as  the  smallest  elements. 

3.1  Example  of  SRI  Alternative 

If,  as  in  the  example  of  Section  2.  6.  2  with  PCFLOS  =  0.  168,  we  want  to  find 
the  probability  of  a  1-km  clear  interval  in  track  of  50  km,  the  SRI  formula  gives 


\ 
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4.  05  km 


M  = 

A  =  20.  10  km 

PI(s'  =  1  km,  T  =  50  km)  =  0. 96  . 

That  is,  the  SRI  formula  estimates  96  percent  probability  of  a  1-km  interval  in 
50  km,  compared  with  the  84  percent  probability  of  the  3D-BSW  model. 

3.2  Comparing  the  SRI  and  3D-BSW  Procedures 

First,  the  number  of  parameters  used  to  characterize  a  cloud  cover  differs 
between  SRI  and  3D-BSW.  In  the  SRI  formula  the  amount  of  cloud  cover  is  essen¬ 
tially  the  only  parameter.  In  the  3D-BSW  model  there  is  the  additional  parameter 
of  scale  distance,  which  is  found  to  vary  with  time  of  day  and  time  of  year,  as  well 
as  geographical  location  and  altitude  above  the  ground.  Eventually  scale  distance 
may  be  found  to  vary  with  cloud  types,  although  this  has  not  yet  been  studied.  The 
data  used  by  SRI  are  best  fitted  by  the  3D-BSW  model  when  the  scale  distance  is 
chosen  to  be  0.  109  km,  or  wavelength  28  km.  In  previous  work  on  clouds,  rainfall 
and  radar  echoes,  we  have  found  scale  distances  of  an  order  of  magnitude  greater. 
This  has  prompted  a  further  comparison  of  results  of  the  SRI  model  with  the  3D-BSW 
model  (Table  5).  If  the  nature  of  the  clouds,  especially  in  horizontal  persistence, 
as  represented  by  the  scale  distance  (r),  should  differ  from  the  clouds  sampled  by 
SRI,  then  the  estimated  likelihood  of  seeing  through  breaks  in  the  clouds  could  differ 
dramatically. 

Second,  the  algorithmic  solution  for  the  SRI  model  is  simpler  than  for  the 
3D-RSVV  model.  However,  once  the  software  program  is  assembled  for  the  latter 
procedure,  it  should  make  little  difference  to  the  computer. 

Third,  in  deriving  their  equations,  the  SRI  authors  assumed  independence  of 
cloud  cover  between  successive  non -overlapping  segments  of  the  track.  In  the 
3D-BSW  model,  on  the  other  hand,  a  field  of  correlation  is  deliberately  constructed 
into  the  model,  with  the  consequent  effect  of  reducing  the  likelihood  of  a  clear  run 
when  the  mean  cloud  cover  is  large  (for  example,  9/10).  This  effect  is  noticeable 
in  Table  5. 
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Table  5.  The  Estimates  of  Probability  of  a  Clear  Interval  of  Length  s'  (km) 
Over  a  Track  of  50  km.  When  the  Sky  Cover  is  0.  9,  for  Which  CFLOS 
Overhead  Probability  is  0.  168  (scale  distance  (r)  is  varied  from  0.  05  km 
to  5  km) 


Interval 

3D-BSW  Model 

Length 
s'  (km) 

SRI 

Model 

r  =  0.  05  km 

r  =  0.  11  km 

r  =  1  km 

r  =  5  km 

0.  5 

0.  988 

0.  977 

0.  883 

0.  35 

0.  32 

1 

0.  96 

0.  93 

0.  83 

0.  34 

0.  27 

5 

0.  51 

0.  26 

0.  45 

0.  32 

0.  21 

10 

0.  17 

0.016 

0.  16 

0.  25 

0.  21 

50  km 

negligible 

negligible 

negligible 

0.  05 

0.  13 

4.  PARAMETER  DETERMINATION 


The  3D-BSW  model  depends  basically  on  two  parameters.  In  this  sense  it 
resembles  most  statistical  models,  in  which  the  most  likely  parameters  are  mean 
and  standard  deviation.  In  the  3D-BSW  model,  the  mean  is  used  but  standard 
deviation  is  replaced  by  the  scale  distance  (r)  or  its  alternate  the  wavelength  (A). 

While  the  mean  or  single-point  probability  is  relatively  easy  to  obtain,  the 
scale  distance  usually  lakes  considerable  effort  that  must  make  the  most  of 
climatological  records.  Accordingly,  an  in-depth  description  of  how  these  param¬ 
eters  are  obtained  with  the  3D-BSW  model,  is  postponed  to  a  later  report. 

A  description  of  parameter  determinations  can  be  found  using  the  2D-BSW 
4  9 

model.  Gringorten's  report  on  cloud  distributions  in  the  vertical  finds  scale 
distance  at  Bedford,  MA  sky  cover,  ranging  from  0.  5  km  at  noon  in  summer 
to  10  km  at  midnight  in  winter.  World-wide,  Burger's1  atlas  of  sky  cover  presents 
similar  values,  which  show  considerable  variation  with  geography.  An  earlier 
paper  showed  that  scale  distance  is  of  the  order  of  1  to  5  km  for  radar  echoes, 

4  to  10  km  for  24-hr  rainfall  in  New  England. 


9.  Gringorten,  1. 1.  (1982)  Climatic  Probabilities  of  the  V  ertical  Distribution  of 
Cloud  Cover.  AFGL-TR-82-0078,  ADA  1  18753. 

10.  Burger,  C.F.  (1985)  World  Atlas  of  Total  Sky  Cover.  AFGL-TR-85-0198, 

ADA  170474. 
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5.  CORRELATIONS 


As  described  above,  the  BSW  model  was  developed  to  provide  useful  simulations 

of  the- weather,  particularly  clouds  and  to  provide  the  probabilities  of  events  or 

conditions  in  3D-space.  The  field  of  correlation  coefficients  (cc)  associated  with 

the  model,  while  not  a  goal  in  itself,  is  an  important  indicator  of  the  validity  of 

the  model,  and  may  be  used  as  a  tool  in  achieving  the  primary  goals. 

While  studies  of  the  spatial  field  of  correlation  have  not  been  as  numerous  as 

studies  of  time-lapse  correlation,  they  have  not  been  neglected.  Buell^1  found 

patterns  in  the  cc's  of  horizontal  winds  aloft,  as  functions  of  distance  between 

12 

stations.  Bertoni  and  Lund  published  some  revealing  diagrams  of  winter  cc's, 
in  the  horizontal,  of  pressure,  temperature  and  air  density  from  the  surface  to 
altitudes  as  high  as  16  km.  For  temperature  they  found  horizontal  cc's  to  decrease 
to  zero  at  distances  of  roughly  1,000  nm,  then  become  negative  to  nearly  -0.2, 
then  back  to  zero  at  roughly  2,  200  nm  (Figure  7), 

5.1  Correlation  in  the  3D- BSW  Model 

The  cc  of  the  END'S  at  two  points  is  to  be  found  as  a  function  of  their  separation 
(s).  The  vertical  scale  is  multiplied  (see  above)  by  a  factor  q  (=  50.  say)  to  make 
the  vertical  measure  cause  a  decay  in  correlation  more  rapid  than  the  horizontal 
measures. 

52,  Derivation  of  the  Analytical  Expression 
for  the  cc 

The  END  (y),  at  any  point  (u,  v,  w)  in  a  three-dimensional  space,  is  a  standard¬ 
ized  average  of  a  set  of  N  heights  (x.,  i  =  1,  N)  of  N  wave  formations,  thus: 

y  =  (Sx.  -  N.x)/«T  .  n/N)  (33) 

where  x  is  the  true  mean  of  all  x^,  o ^  is  the  true  standard  deviation  of  x^.  In  the 
3D- BSW  model,  the  x's  have  a  rectangular  distribution  of  values  between  0  and  1, 
hence 

x  =  1/2 

(34) 

a  =  1N12 

X 

11.  Buell,  C.E.  (1962)  Two-Point  Variability  of  Wind,.  Final  Report,  AF19(604)- 

7282,  AFCRL-62-889,  Kaman  xNuclear,  Colorado  Springs,  Colorado. 

12.  Bertoni,  E.  A. ,  and  Lund,  I.A.  (1964)  Winter  Space  Correlations  of  Pressure, 

Temperature  and  Density  to  16  km.  Environment  Res.  Papers  (No.  75), 
AFCRL-64- 1020,  ADA  611002,  29  pp. 
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to  make 


y  =  (Ex.-  N/2)/<Nfi57n\ 


5.5  km  O'STamcE  ihuno»eos  of  nautical  miles i 


Figure  7.  Horizontal  Decay  of  Temperature  Correlation  as  a 
Function  of  Distance  (nm)  Between  Stations  12 


In  any  one  wave  formation  the  cc  between  Xj  and  x 2>  that  are  separated  by  the 
distance  s,  is  given  by 


px(s>  = 


E(x1  x2)-x  2 


(36) 


In  the  set  of  N  wave  formations,  the  cc  between  y^  and  y2«  at  any  two  points 
separated  by  distance  s,  is  given  by 


pj  s>  =  E<y1y2) 


(37) 


A  brief  manipulation  of  these  equations  will  reveal  that 


Py(s>  =  Px<s>  • 


(36) 


Hence,  to  find  p^(s)  it  is  sufficient  to  find  Px(s),  now  simply  denoted  as  p(s). 

It  is  convenient,  at  this  stage,  to  measure  all  distances  in  units  of  the  wave¬ 
length  (A).  That  is,  in  each  wave  formation  the  crests  will  be  one  unit  of  distance 
apart  (A  =  1  temporarily). 

Without  loss  of  generality  we  can  limit  our  attention  to  one  wave  formation  and 
think  of  the  wave  crests  as  thin  flat  sheets  of  plywood  separated  by  one  unit  of  dis¬ 
tance,  all  sheets  parallel  to  the  V-W  plane.  In  Figure  8  the  U,  W-axes  are  in  the 
flat  surface  of  the  paper,  the  V-axis  and  the  plywood  sheets  rise  out  of  the  paper. 

We  place  the  first  point  (C)  along  the  U-axis  at  distance  Xj  from  the  leading 
edge  of  the  middle  wave,  and  restrain  x^  to  vary  uniformly  from  zero  to  1.  0.  In 
the  2D-BSW  model  the  wave  height  at  C  was  its  phase  (x^).  In  the  3D-BSW  model 
the  locus  of  the  second  point  (P)  at  distance  (s)  from  the  point  C  is  a  spherical 
surface  of  radius  (s),  represented  only  by  the  circle  in  Figure  8.  The  height  (x2) 
of  the  second  point  at  P  is  given  by 


x2  =  X1  +  s”  cos  ®  +  ® 


(39) 


where  9  is  the  angle  between  the  line  CP  and  the  U-axis,  and  where 


6  =  -1,  -2, - if  P  is  in  the  first,  second, - right-hand  wave. 


0  if  P  is  in  the  middle  wave  (as  shown  in  Figure  8), 


=  1,2, - if  P  is  in  the  first,  second, - left-hand  wave. 
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Figure  8.  A  Cross-section  of  a  Wave  f  ormation  in  the  U-W  Plane  to 
Illustrate  the  Relation  Between  Two  Points  (C.P)  Separated  by  Distance 
(s).  The  wave  phase  resembles  a  flat  sheet  extending  into  the  V -dimension 
(into  and  out  of  the  paper) 


The  covariance  of  (Xj.  x^)  is  obtained  by  finding  the  average  of  (Xj  x.,)  over 
the  sphere  of  radius  s,  whose  surface  area  is  47Ts^.  The  element  of  area  at  P 
is  given  by 


6A  =  (s,60  ).  (s.  sin  £  .  6  £)  (40) 


where  £  is  an  angular  measure  around  the  small  circle  of  radius  s.  sind  ,  centered 
on  the  U-axis.  For  integration  we  must  have 


0  S  9  5  n 


0  §  £  £2n  . 
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Hence 


71  2  7! 


E  (x,  O  = 


1 


2  . 


^  /  /  /  x  (x  +  s  cos  0  +  6 )  .  s  sin  6  .  dx.  d  6 


4  js  x  —  0  0  =  0  ?=0 


Integrating  with  respect  to  and  with  some  terms  vanishing, 

1  71 

2 

2E  (x.  x„)  =  ,  +  /  xdx  /  6  .  sin  8  dd  . 

1  J  x=  0  0=0 


When  0  s  =  1.  6  is  restricted  to  6  =  -1,0,  1.  Hence 

s  77  1  01 


2E 


(x1x9)  =  —+  /  xdx  /  sin  0 .  d0  +  /  xdx  /  (~1) 

^  ^  A  _  a  «  A  A 


x  =  0 


8  =7!  -8 


x=  1  -  s 


0  =0 


where 


cos  0  j  =  ( 1  -  x)  /s 


cos  (t;  -  0  2)  =  x/s 


which  simplifies  to 


1 2.  E  (Xj  x^)  =  4  -  3s  +  2s 


whence,  by  Eqs.  (34)  to  (36) 
p(s)  =  1  -  3s  +  2s2 


for  0  i  s  %  1 


When  1  <  s  S  2,  then  6  can  have  two  or  more  values.  Or 
6  -  -2,  -1,  0,  1  or  2  . 

The  solution  for  E(Xj  x^)  is  much  longer  and  involved,  but  equally  valid, 
finally 


p(s)  =  (1  -  3s  +  2s2)  -  6  (s  -  l)2/s  for  1  <  s  =  2  . 


.d  C  . 
(41) 


(42) 


sin  0  d0 


(43) 


(44) 


It  gives. 


(45) 
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Figure  9  is  a  plot  of  p(s)  vs  s,  where  s  is  in  units  of  the  wavelength.  The  cc 
becomes  zero  at  one-half  wavelength,  becomes  negative  with  minimum  -0.  12, 
then  zero  again  at  1,  1.5  and  2.0  wavelengths.  As  s  increases,  the  value  of  pis) 
will  oscillate  from  positive  to  negative,  and  will  become  vanishingly  small. 

When  A=  300  km,  the  cc  very  nearly  equals  0.  T9  over  the  unit  distance  of 
1  km.  For  other  A  we  make  cc  a  function  of  distance,  (s'  km),  thus: 

p(s')  =  1  -  3s  +  2s^  where  s  =  s' /A  for  0  £  s'  =  A  .  (46) 

The  value  of  A  that  will  minimize  the  root  mean  square  difference  between  the  cc's 
of  the  previous  Model  B  and  the  3D-BSW  Model,  for  distances  up  to  the  distance 
at  which  Model  B  yields  cc  equal  to  zero,  is  close  to  260  km.  We  choose  to  set 
A^  =  2^  =  256  km.  The  cc  will  drop  to  zero  at  s'  =  2^  =  128  km;  the  cc  at 
s'  =  1  km  will  be  0.  9883 - . 

5.3  Correlation  in  the  Vertical 

The  climatic  records  (WBAN10  and  TDF14)  give  information  on  the  amount  of 

cover  (in  tenths)  at  one  to  four  levels,  plus  the  height  of  the  cloud  bases.  But  how 

1 3 

can  we  derive  cc's  from  this  information?  The  Air  Weather  Service  nephanalvsis 
model  3-D  Neph,  gives  3-hourly  information  on  cloud  presence  at  some  15  levels 
in  the  atmosphere,  from  which  tables  of  the  frequency  distribution  of  cloud  cover, 
in  layers,  have  been  compiled.  The  record  might  be  explored  to  see  if  cc's  also 
could  be  computed  or  estimated.  We  are  unaware  of  any  published  results,  and 
for  the  present,  we  consider  cc  in  the  vertical  to  be  unknown.  (In  the  earlier 
sections  of  this  report,  the  value  of  q=  50  is  only  a  guess.  ) 

The  application  of  the  3-D  model  to  produce  clouds  in  three  dimensions  should 
bypass  the  cc.  The  approach  to  the  distribution  of  cloud  cover  in  the  vertical  is 
to  treat  a  cloud  layer,  of  thickness  (h),  as  though  it  were  an  entity  with  one  mean 
cloud  cover  (PQ)  and  scale  distance  (r)  for  that  layer.  The  mean  cloud  cover,  in 
layers,  is  obtainable  from  the  3-D  Neph  data.  It  is  also  obtainable  from  the  TDK  14 
information  on  cloud  cover  from  the  ground  to  height  (H)  and  by  using  the  parameter 
G  (h,  H),  previously  described.  * 

For  the  present,  the  scale  distance  of  a  layer  that  has  specific  lower  and 
upper  levels  is  taken  to  be  the  same  scale  distance  as  that  determined  for  the  whole 
layer  from  the  ground  up  to  the  top  of  the  cloud  layer. 


13.  Fye,  F.K.,  Major  USAF,  (1978)  The  AFGWC  Automated  Cloud  Analysis 
Model.  AFGWC  Technical  Memorandum  78-002,  ADA  057176, 

Hq  AFGWC,  Qffutt  AFB,  Nebraska. 
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Relation  of  Correlation  Coefficient  p(s)  to  Distance  (s)  in  the  3D-BSW  Model.  In  this  diagram 


5.4  Probability  of  Joint  Events  at  Two  or 
more  Stations 

The  main  thrust  of  this  report  is  aimed  at  spatial  events  within  a  known  space, 

rather  than  specific  points  within  the  space.  However,  since  the  3D-BSW  model 

provides  the  correlation  coefficient  between  points  of  known  distance  apart,  it  is 

possible  to  obtain  the  probability  of  joint  events,  through  application  of  the  multi- 

14 

variate  normal  distribution. 

The  joint  probability  of  two  points  exceeding  a  threshold  is  given  by  the  integra¬ 
tion  of  the  bivariate  normal  distribution.  Of  the  many  ways  to  calculate  this  integral, 
we  prefer  to  use  Eq.  (31)  of  Owen1'’  along  with  Algorithm  No.  15  of  Yamauti.  16 

For  the  joint  probability  of  three  or  more  points,  algorithms  become  increas¬ 
ingly  complex.  We  find  it  simpler  to  cumulate  the  results  using  a  multitude  of 
simulations. 

6.  ALTERNATIVES  TO  THE  3D-BSW  MODEL 


The  3D-BSW  model  is  basic  to  the  methodology  of  this  report.  It  succeeds 
the  previous  Gringorten  model  B  and  the  2D-BSW  model,  and  is  subject  to  change 
or  improvement  in  the  future.  A  few  alternatives  have  been  considered,  including 
multi-dimensional  models  to  include  characteristics  other  than  the  spatial,  such 
as  time.  Instead  of  sawtooth  waves,  one  model  has  wedge-shaped  waves  and 
another  sinusoidal  waves. 

6.1  The  4D-BSW  Alternative 

If  simulation  is  to  be  realized  by  the  3D-BSW  model  in  space,  coupled  with  a 
process  in  time,  the  procedure  will  call  for  the  stochastic  production  of  a  three- 
dimensional  field  of  END-values,  yt  (u,  v,  w),  at  each  point  (u,  v,  w)  for  an  initial 
time  (t).  Then  another  field  of  END-values  must  be  generated  stochastically,  and 
linked  to  the  initial  field,  to  give  time- related  END-values.  For  every  additional 
time  another  field  would  be  generated,  with  its  attendant  random  numbers.  For  a 
20-year  run  with  a  1-min  time  step,  this  would  require  the  random  generation  of 
over  10  million  fields. 


14.  Gringorten,  1.  1.  (1978)  Conditional  Joint  Probabilities,  AFGI.-TK-78-0238, 

ADA  063817. 

15.  Owen,  D.  B.  (1980)  A  table  of  normal  integrals,  C'ommun.  Statist,  -Simul. 

Computa. ,  J39(4):389-419. 

16.  Yamauti,  Ziro,  Ed.  (1972)  Statistical  Tables  and  Formulas  With  Computer 

Applications,  Japanese  Standards  Association. 


To  avoid  the  dilemma  just  described,  the  3D-BSW  model  has  been  advanced  to 
the  4D-BSW  model,  in  which  the  fourth  dimension  is  devoted  to  time.  The  repre¬ 
sentation  of  time,  however,  is  just  one  of  many  possibilities  for  a  4D-model.  As 
of  this  writing,  the  simulation,  by  the  4D-BSW  model,  of  a  field  of  values  in  one, 
two,  three  or  four  dimensions  has  been  tried  successfully.  The  decay  of  correla¬ 
tion  with  distance  for  one  wavelength  (A)  is  given  by 

p(s)  =  1  -  8s/7i  +  3s^/2  for  s  =  s'/A  %  i  .  (47) 

The  application  of  the  climatic  probabilities  has  been  explored.  The  areal  coverage 
probabilities  have  been  found  to  be  similar  to  the  3D  coverage.  Much  more,  how¬ 
ever,  needs  to  be  done. 

6.2  Sinusoidal-Wave  Model 

The  sinusoidal-wave  model  has  its  attractions  and  may  yet  be  developed 
intensively.  The  primary  drawback  is  the  time  required  to  produce  synoptic  maps. 
Wave  formations  are  produced  stochastically,  as  in  the  3D-BSW  model,  but,  instead 
of  Eq.  (14)  for  the  wave  height  (x)  at  a  point  (u,  v,  w),  the  equation  becomes 

x  =  sin  j 2ir(h  +  D/A)  j.  48) 


it  differs,  also,  from  the  3D-BSW  model  in  the  shape  of  the  curve  of  correlation 
vs  distance  (figure  10).  I'or  s  near  zero,  the  3D-BSW  correlation  decreases  nearly 
exponentially ,  while  the  sine  correlation  has  an  exponential-squared  decay. 

Instead  of  Eqs.  (44)  or  (45),  the  formula  for  cc  vs  distance  becomes 


p(s) 


sin  (27?s) 
2tjs 


for  all  s  g  0  . 


(49) 


While  the  exponential -type  decay  of  cc  serves  us  well  for  meso-scale  phenomena, 
such  as  cloud  and  precipitation  patterns,  the  exponential-squared  type,  which 
f  igure  10  suggests,  may  be  more  appropriate  for  macro-scale  phenomena.  There 
was  a  hint  of  this  preference  when  the  earlier  Gringorten  Models  A  and  B  were 
studied,  in  particular  for  upper-air  temperatures. 
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7.  SUMMARY  AND  CONCLUSIONS 


The  3D-BSW  model  was  developed  to  estimate  probabilities  of  events  (X)  in 
one-,  two-  or  three-dimensional  space.  It  requires  a  parameter  (A)  such  that  the 
correlation  coefficient  of  the  END'S  (y)  of  X  at  two  stations  will  reduce  to  zero  at 
a  separation  of  (A/ 2)  km  between  the  stations. 

Symbolically  the  algorithms  have  been  written  to  give: 

(a)  PA  (X,  A,  F,  r),  the  probability  that  a  threshold  condition  (X) 
will  be  exceeded  only  in  a  fraction  (F/10)  or  less  of  the  area  (A), 

(b)  PL  (X,  s',  F,  r),  the  probability  that  X  will  be  exceeded  only  over 
(F/10)'s  or  less  of  the  line  of  length  (s'), 

(c)  PI<X,  s',  T,  r>,  the  probability  that  there  is  a  distance  (s')  out 
of  an  overall  line  of  travel  <T),  over  which  the  condition  (X)  is 
not  exceeded  (that  is,  the  minimized  maximum). 

The  3D-BSW  model  has  special  application  to  cloud  cover.  If  the  threshold 
is  considered  to  be  the  no-cloud  condition,  then  the  algorithms  provide  the  prob¬ 
abilities  of  a  clear  condition  in  the  fraction  (F)  of  an  area  (A)  or  line  of  length  (s'), 
along  with  the  probability  of  a  cloud-free  interval  of  length  (s’)  in  a  total  line  of 
travel  (T). 

Methods  for  the  determination  of  the  parameter  (A)  or  (r)  are  still  to  be 
published.  The  data  sources  are  primarily  the  climatological  summaries  of  hourly 
data. 

As  might  be  expected,  the  fields  generated  bv  the  model  are  characterized  by 
correlation,  varying  with  the  distance  of  separation  between  points.  A  vertical 
distance,  with  respect  to  correlation,  is  comparable  to  a  much  greater  distance  in 
the  horizontal.  This  relationship  needs  to  be  studied  further. 

Paradoxically  the  determination  of  the  joint  probability  of  a  threshold  condition 
at  three  or  more  points  is  more  difficult  than  finding  such  probability  for  a  whole 
a  re  a. 

The  3D-BSVV  model  may  be  improved,  or  replaced,  eventually.  So  far,  it  has 
emerged  as  the  optimum  tool  for  simulating  synoptic  images  and  climatic  frequencies 
of  conditions  or  events  in  space. 
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Purpose: 


Comment: 


Step  0. 


Appendix  A 

Siwtatnf  a  1-  m  2-OmmmmmI  FM  af  END'S 


To  produce  stochastically  a  one-  or  two-dimensional  field  or  portrayal 
of  END -values  of  a  weather  element  using  the  3D-BSW  model.  A 
one -dimensional  field  is  produced  by  keeping  two  components 
(u,  v  or  w)  constant,  a  two-dimensional  field  is  produced  by  keeping 
one  component  constant. 

Data  are  not  required.  Instead  a  software  program,  to  generate 
random  numbers,  each  from  a  rectangular  distribution  of  numbers 
from  0  to  1,  is  used.  Each  map  or  configuration  will  require  the 
generation  of  several  dozen  random  numbers. 

Choose  A;  a  wavelength,  the  parameter  of  the  3D-BSW  model. 

Choose  R  :  an  initial  random  number  or  seed.  Note:  It  has  been 

o  - 

found  helpful  to  record  this  number,  since  it  may  be 
desired  to  intensively  explore  a  small  part  of  the 
synoptic  field  in  a  supplementary  investigation. 

Choose  L:  the  number  of  wave  formations  per  map. 

Choose  (u  ,  v  ,  w  )  the  starting  point. 

Choose  Mu:  the  number  of  intervals  to  the  side  in  the  u-direction 

Choose  M  :  the  number  of  intervals  in  the  v-direction. 

v  • 

Choose  M  :  the  number  of  levels  in  the  w -direction, 
w 

Choose  6s1  (km):  the  interval  in  the  horizontal. 

Choose  6w'  (km):  the  vertical  distance  between  levels. 

Choose  q:  the  quotient  to  "expand"  the  vertical  scale. 
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I 


« 


Note:  The  following  description  of  the  method  is  strongly  influenced  by  the  BASIC 
language. 

Step  1.  For  I  ~  1  to  L  (for  generating  the  Ith  wave  formation) 

Find  R^,  a  random  number  between  0.  0  and  1.  0 

Set  h„  =  R, 
f  h 


Step  la.  To  find  cos  a cos  f}^,  cos  y^. 


Find:  R  ,  R  ,  R  ,  three  random  numbers  between  0.  0  and  1.  0 
u  v  w 


Find:  r  =  R 
u  u 

r  =  R 
v  v 

r  =  R 
w  w 

Find  (OP)2  =  r2  +  r2  +  r2. 

u  v  w 


-  1/2. 
-  1/2. 
-  1/2. 


If  (OP)  >  1/4,  then  reject  the  three  random  numbers,  and  repeat  Step  la. 
Else,  when  (OP)2  %  1/4,  find  cos  =  ru/(OP) 

cos  /3g  =  rv/(OP) 


cos  7 


rw/(OP>  . 


Continue  Step  1,  as  long  as  (  <  L. 

Step  2.  For  m  =  0  to  M  (for  the  m  th  level) 

r  WWW 

Set  w  =  (w  +  m  *  6w’)-  q. 
o  w 

Step  3.  For  my  *  0  to  (for  the  mvth  point  in  the  V -direction). 

Set  v  =  v  +  m  .6s'. 
o  v 

Step  4.  For  mu  =  0  to  Mu  (for  the  m^th  point  in  the  U -direction). 

Set  u  =  u  +  m  ■  6s’. 
o  u 

Initialize  Z  (u,  v,  w)  =  0. 

Step  5.  For  f  =  1  to  L  (for  the  {th  wave  formation). 

Find  D  =  u.  cos  +  v  •  cos  /Sj  +  w  cos  >  ^. 

Find  T  -  h^  +■  D/A  +  I  (where  I  is  a  sufficiently  large  integer, 

to  make  T  positive;  if  BASIC  is  used,  I  isn't  necessary). 

Find  (u,  v,  w)  =  FRA  (T >  =  T  -  INT(T). 

Note:  Section  2.  1,  5  describes  how  to  find  the  wave  height  (x),  tailored  to  both 
FORTRAN  and  BASIC. 
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L 

Au  i  x.tu.  v,  w)  to  J2  x.  (u,  v,  w) . 
f  f  =  1 

Continue  Step  5  as  long  as  I  f  L. 

L 

Find  y(u,  v,  w)  =  V ( 12L) •  [  -r-  S  x.(u,  v,  w)  -  1/2]  . 

L  f  =  1  * 

Continue  Step  4,  as  long  as 

Continue  Step  3,  as  long  as  mv  =  Mv- 

Continue  Step  2,  as  long  as  m  %  M  . 

r  b  w  w 

(This  exercise  provides  END-values  at  all  specified  grid  points,  in  space.  Since 
this  appendix  was  written,  a  modification  of  the  3D-BSW  model  has  enabled  us  to 
use  one  of  its  three  dimensions  to  closely  resemble  the  O-U  process  for  time, 
which  will  be  described  in  a  later  report.  ) 


Appendix  B 


Estimating  ths  Probability  That  a  Threshold  Valaa  WiN  bt 
Eicatdad  in  Only  (F/10)ths  of  a  Sqoara  Araa  (A) 


Purpose:  To  estimate  PA  (y  ,  A,  F.  r),  the  probability  that  a  threshold  value 

of  X,  or  its  END  (y,  ),  will  not  be  exceeded,  except  in  (F/10)ths  or 
less  of  the  square  area  (A);  scale  distance  is  (r). 

Algorithm: 

Step  0:  Begin  with:  threshold  (y  )  or  the  single-point  cumulative 

probability  (P  )  o"  X,  plus  scale  distance  (r  km),  plus  area 
°  2 

of  floor  space  (A  km  '). 

Step  1.  Compute  entries  yQ  (if  necessary)  and  z. 

Find  z  =  (in  v/A/r)/ln  2. 

If  given  Pq,  find  y^.  The  recommended  formula  is  the  NBS  aoproximation; 

yo  =  k  [t  -  (aQ  <-  aj.  D/d  +  bjt  +  bgt2)]  (Bl) 

where 

aQ  =  2.307  53 
a:  =  0.  27061 
bj  =  0. 99229 


0. 04481 


and  where 


k  =  -1,  t  =  \/(fn  1/p2)  for  p  =  Pq  s  1/2 

k  =  1,  t  =  *J(fn  1/(1  -  p)2)  for  p  =  Pq  >  1/2  . 

Step  2.  To  solve  for  the  probability  of  the  fractional  cover  (F/10): 

Note;  In  this  algorithmic  solution,  F  need  no  longer  be  an  integer. 

If  0  %  F  %  5,  substitute  F'  *  F,  y  1  =  y  . 

If  5  <  F  =  10,  substitute  F*  =  10 -F,  y  1  =  -y  . 

The  END,  y(F'),  is  obtained  by  interpolation  between  y  jlNT(F')} 
and  y  jlNT(F')+  1  }  .  Both  of  these  latter  values  are  obtained  by 

y(F')  =  a  +  0  i-  y  t  5  (B2) 

where  Ct,  fi,  y,  5  are  solved  by  the  equations  in  Table  Bl,  supported 
supported  by  Tables  B2  and  B3. 

The  interpolation  procedure  is  aceomolished  by 

y(F')  =  y|lNT(F')(  f  |f'-INT(F')|  j  y[INT(F')+  1] -y  [INT(F«  )]| .  (B3> 

Find  the  probability,  P  corresponding  to  the  End,  y(F'). 

The  recommended  formula  is  the  NBS  approximation 

P  =  f +m  [2(1+c1x+c2x2+c3x3+c4\4)4]'1  (B4> 


where 

Cj  =  0.  196  854 

c2  =  0. 115  194 

c3  *  0.000  344 

c.  =  0. 019  527 
4 

x*  ly(F')| 

1  =  0,  m  =  1  for  y(F')  %  0 
f  *  1,  m  =  -1  for  y(F' )  >0. 

For  0  I  F  1  5,  set  PA  (yQ,  A,  F,  r)=  P 
For  5  <  F  %  10,  set  PA  (y  ,  A,  F,  r)  =  1  -  P. 
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should  be  read  as  y  P  as 


Table  B2.  Constants  for  the  Expression 
(0)  in  Table  B1 


n 

a 

n 

b 

n 

C 

n 

0 

0. 0900 

0. 0000 

0.00 

1 

-0.  007  1 

0.0086 

0.  53 

2 

-0.  0205 

0. 0097 

0.  52 

3 

0. 0229 

0.  0045 

0.  64 

4 

-0.  0517 

0. 0224 

0.  45 

5 

0. 0260 

0.  0087 

0.71 

Table  B3.  Constants  for  the  Expression  (6)  in  Table  B1 


m  =  1 

....  .  . 

m  =  2 

F 

d<F)m 

f(F)™ 

m 

g(F) 

&  m 

d(F) 

m 

f(F)m 

g(F) 

6  m 

0 

0.  1420 

-2.  4090 

-0.  44 

0.  1908 

-3.  8406 

-0.  56 

1 

1.  0686 

-2.  2854 

-0.  46 

1.  87  10 

-4.  2412 

-0.  55 

2 

1.  9440 

-2.  7842 

-0.  41 

3. 0540 

-4. 4752 

*0.  39 

3 

2.  2560 

-2.  7356 

-0.  40 

3. 3960 

-4. 5228 

1 

© 

09 

o 

* 

2. 0328 

-2.  2910 

-0.  41 

2. 0390 

-2. 5702 

-0.  7  5 

5 

1.  2356 

-1.  2300 

-1.  30 

1.  6124 

-1.  5434 

-1.  51 

The  following  BASIC  program  (filename  B:PYOA  FR.  BAS)  estimates  PA  (y  ,  A,  F,  r  ). 
It  requires  the  inputs: 

scale  distance,  entered  as  R  (km), 

2 

area  of  floor  space,  entered  as  A  (km  ), 
fraction  of  area,  entered  as  F, 

the  cumulative  probability  of  the  threshold  (X)  entered  as  P  , 
and  yields  the  output  (P),  the  probability  of  the  event. 
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10  REM  B» PYOAFS.  BAS  ON  GRINGOKTEN  DISK  NO.  4.  IB  DEC  1986 

20  SEN  TO  FIND  THE  PROBABILITY  THAT  X.  ITS  END  <Y0>  WHOSE  CLIMATIC  FREQUENCY 
IS  PO.  HILL  BE  EXCEEDED  IN  !F/10.)THS  OR  LESS  OF  AREA  A:  PARAMETER  IS  SCALE 
DISTANCE  (R* 

30  REtf  PO  HEREIN  HILL  BE  CUMULATIVE  PROBABILITY 

40  REM  THE  FOLLOWING  ARE  CONSTANTS  FOR  THE  ALGORITHMS 

50  AO  *  2.30753 

60  A1  =  .  27061 

70  B1  «  .99229 

80  B2  *  .04481 

90  Cl  =  .196854 

100  C2  =  .115194 

110  C3  *  .000344 

120  C4  «  .019527 

121  OPTION  BASEO 

130  DIM  AN<6>.  BN (6 ) .  CN<6! 

140  DIM  DFM<  2.6) 

150  DIM  FFM<2.6> 

160  DIM  GFM (2.61 
170  DIM  RD ( 2.  2 ) 

180  FOR  I  *  6  TO  5 

190  READ  AN ( I ) .  BN < I ) .  CN(I? 

200  NEXT  I 

210  DATA  0,  0.  0. -.  0071.  .0086.  .53.  -.0205.  .0097.  .52 

220  DATA  .0229.  .0045.  .64.  -.0517,  .0224.  .45,  .0260,  .0087,  .71 

230  FOR  H  =  1  TO  2 

240  FOR  J  =  0  TO  5 

250  READ  DFH(M.J).  FFM(H.J).  GFM'H.J) 

260  NEXT  J 
270  NEXT  M 

280  DATA  .1420.  -2.4090.  -.44.  1.0686.  -2.2854,  -.46 

290  DATA  1.9440.  -2.7842.  -.41.  2.2560.  -2.7356,  -.40 

300  DATA  2.0328.  -2.2910.  -.41.  1.2356.  -1.2300.  -1.30 

310  DATA  .1908.  -3.8406.  -.56.  1.8710,  -4.2412,  -.55 

320  DATA  3.0540.  -4.4752.  -.39.  3.3960.  -4.5228.  -.50 

330  DATA  2.0690.  -2.5702.  -.75,  1.6124.  -1.5434.  -1.51 

340  REM  STEP  0  FOR  THE  INPUTS 

350  INPUT  "SCALE  DISTANCE  =":R 

360  INPUT  "AREA*”: A 

370  INPUT  "FRACTION  =  “ :  F 

380  INPUT  -P0=“iP0 

390  REM  STEP  1  TO  TRANSFORM  THE  VARIAELES  INTO  ALGORITHM  ENTRIES 

400  REM  TO  CALCULATE  2 

410  TEMP  =  SCR ( A '  R 

420  2  =  LOG (TEMP  > • L0G<  2 > 

430  REM  TO  CALCULATE  YO  FROM  PO 
440  IF  POO. 5  THEN  480 
450  K  =  -1 

460  T  =  SCR ; LOG  < 1  - P 0  *•> 

470  GOTO  500 

4  80  K  =  1 

490  T  =  SOR ( LOG  1 1 /  1 - P 0 ;  2, 

500  YO  =  K»(T  -  'AO  ♦  A1*T*  •’  1*B1*T  ♦  E2*T  2  i 
510  REM  TC  SUBSTITUTE  Y  0  F  =  -  Y  :  AND  FF  =  10-F  WHEN  F  20  .5 
520  IF  F‘  5  THEN  560 
530  FP  =  F 

5  4  0  Y  0  P  =  Y  0 
550  GOTO  590 
560  FP  =  10-F 
570  YOF  =  -YO 
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580  REM  STEP  2  TO  CALCULATE  THE  END  <Y)  OF  THE  DESIRED  PEOEAEILITY 

590  IF  ::1  THEN  650 

600  REM  FOR  THE  CASES  WHEN  Z'l 

610  ALPHA  =  Y OP  +  ( . 006»EXP !  .  8*Z ) -  .  003 1 » < FP- 5 • 

620  YP  =  ALPHA 
630  GOTO  1410 

640  REM  FOR  THE  CASES  WHEN  Z=:i 
650  IF  207  THEN  900 

660  REM  FOR  THE  CASES  WHEN  1=CZ=07:  FIRST  WE  FIND  BETA 
670  FIRST  =  (.9981  ♦  .  00 11*  EXP  <  .£■-'<  Z  '  !  •»  YOF 
680  N=  5  -  INTIFP) 

700  TEMP  =  EXP ( CN ( N ) *2 ) 

710  SECOND  =  (AN(N)*BN<N)»TEMP'*<FP+N-6> 

720  SUM  =  0 

730  IF  N  =  0  THEN  790 

740  NMIN1  =  N  -  1 

750  FOR  I  =  0  TO  MM INI 

770  SUM  =  AN ( I '  +  EN (I ) »EXP (CN( I » *Z *  *  SUM 
-’8  0  NEXT  I 

790  EETA  =  FIRST  *  SECOND  -  SUM 
900  REM  TO  fIND  THE  GAMMA  TERM 
810  IF  Z  =05.5  THEN  880 
9  20  IF  FP  ;•  4  THEN  990 
830  IF  YOF  !2.5  THEN  830 

840  REM  FOR  THE  CASE  WHEM  5  .Z  0  ?  >.  YOF.  2.  5 

850  GAMMA  =  (  .  0  01"7*EXP  !  .  95  «  Z  >  -  .  3129  i  «  <  .  -  .  7  *  Y  OP  *  1  1  -  .  25*FF 

860  YP  =  EETA  *  GAMMA 

870  GOTO  1410 

980  YP  =  EETA 

890  GOTO  1410 

900  REM  FOP.  THE  CASE  WHEN  Z!  7 
910  REM  WE  FIRST  FIND  Y!Z=7) 

920  FIRST  =  ( .9981*.  0'.11*EXP'  6.  23  ■  1  *Y0P 
930  N  =  5  -  I  NT ( FP 1 
950  TEMP  =  EXP'CN<N)»7. 

960  SECOND=  t  AN ( N ;  ♦  BN i N ) « TEMP > » ( FP * N - 6 i 

970  SUM  =  0 

980  IF  N=0  THEN  1040 

990  NMIN1  =  N - 1 

1000  FOR  I  =  0  TO  NMIN1 

1020  SUM=  SUM  +  AN ( I )  *  BN ( I ) »EXP < CN i I > *" > 

1030  NEXT  I 

1040  BETA  =  FIRST  +SECOND  -  SUM 
1050  IF  FP >4  THEN  1100 
1060  IF  Y0PC2.5  THEN  1100 

1070  GAMMA  =  (  .  0017*EXP  ;6.  65' -  .  31 29  >  *  .  .  -’5- .  3*  VC  P  '•<  C 25  *FP 

1080  YP7  =  BETA  ♦  GAMMA 
1090  GOTO  1110 
1100  YP7  =  BETA 

1110  REM  FOR  THE  ADEITONAL  TERMS  UHEN  7 
1120  IF  FP=5  THEN  FP=4. 999999 

1130  REM  FP  SHOULD  EE  JUST  LOW  ENOUGH  TO  AVOID  ROUNDOFF  TO  5  EY  THE  COMPUTER 
1140  M  =  INT (Z-5. 5) 

1150  REM  M  IS  EITHER  1  OR  2 
1160  RD! 1. 1  )  =  2-7 
1170  RD ( 2. 1 )  =  .5 
1180  RD  (  2.  2  )  =  Z--\5 

1190  REM  EEC AUSE  F  MAY  NOT  DE  A  WHOLE  NUMBER  THE  FOLLOWING  PROCEDURE  IS  USEE 
120'  FLOWER  =  INTIFP > 

1210  FUTFER  =  FLOWER  *  1 
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1220  FOR  FI  =  FLOWER  TO  FUPPER 

1240  REM  FOR  THE  CORRECTIVE  TERM  WHEN  Z>7 

1250  TERM  =  YP7 

1260  IF  M=2  THEN  1310 

1270  REM  THE  CASE  WHEN  M=1 

1280  TEMP  =  EXP  < GFM ( 1, FI >  *Y0F ) 

1290  TERM  =  TERM  +  ( DFM ( 1 .  FI > +FFM ( 1 , FI ) *TEMP ) *RD ( 1 , 1 ) 

1300  GOTO  1360 

1310  REM  THE  CASE  WHEN  M  =  2 

1320  FOR  J  =  1  TO  2 

1330  TEMP  =  EXP ( GFM  <  J, FI >  * YOP  > 

1340  TERM  =  TERM  +  ( DFM < J,  FI > +FFM ( J,  FI ) *TEMP ) *RD ( 2, J > 

1350  NEXT  J 

1360  IF  FI  =  FUPPER  THEN  1390 
1370  YFIRST  =  TERM 
1380  NEXT  FI 
1390  YSECOND  =  TERM 

1400  YP  =  YFIRST  +  <  FF-FLOWER) * t YSECOND  -  YFIRST) 

1410  REM  CONTINUE 
1420  IF  FI 5  THEN  Y  =  -YP 
1430  IF  F  = 1 5  THEN  Y  =  YP 

1440  REM  STEF  3  TO  CALCULATE  THE  DESIRED  PRODAE ILITY 

1450  IF  YIO  THEN  1500 

1460  L  =  0 

1470  M  =  1 

1480  X= - Y 

1490  GOTO  1530 

1500  L  =  1 

1510  M  =  -1 

1520  X  =Y 

1530  DEN  =  ( 1  +  C 1  * X  +  C 2 * X  2  +  C 3 * X  3  +  C 4 « X  '  4  )  4 
1540  P  =  L  +  M/ ( 2* DEN  > 

1550  PRINT  "FOR  FRACTIONAL  COVER  F=":F;’\  CLIMATIC  PROBABILITY  P0=” ; PO; " IN  AR] 
A= “ ; A,  "  P=  “  ; P 
1560  GOTO  380 
1570  END 
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Purpose: 


Algorithm: 
Step  0. 


Step  1. 


Step  2. 


Appendix  C 

Estimating  tha  Probability  That  a  TbraahoM  Vakia  (y,)  Will  ba 
Excaadad  ia  Only  (F/10)tha  of  a  Liaa  Langtfc  (s') 


To  estimate  PLfy^.  s’,  F,  r),  the  probability  that  the  threshold  of  X, 

or  its  END  (y  ),  will  not  be  exceeded,  except  in  (F/  10)ths  or  less 
-  o 

of  the  line  of  length  (s');  scale  distance  is  given  (r). 

Begin  with: 

Single-point  cumulative  probability  (Pq)  or  its  END  (y^), 
scale  distance  (r  km)  or  wavelength  (A  =  256  r  km),  line 
length  (s'  km). 

Compute  the  entries:  y  (if  necessary),  and  z. 

Find  z  =  In(s'/r)/ln  2.  (Cl) 

Find  yQ;  the  recommended  procedure  is  in  Appendix  B. 

To  solve  for  the  probability  of  the  fractional  cover  (F/10)ths  of 
the  length  (s'):  Note:  F  is  not  necessarily  an  integer. 

For  0  §  F  =  5,  substitute  F'  =  F,  y,  '  =  y  • 

For  5  <  F  5  10,  substitute  F'  =  10-F,  y  '  =  -v  . 

’  o  o 

Fnd  T?r  nr 

where 

■  An  *  B„V  '  S/o'2  K2) 
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where 


An  =  aa+bn‘  exP(cn2) 

Bn  =  da+fnexp(Vz) 

C  =  h  +  q  .  exp  (r  .  z)  for  z  >  6 
n  nnn 

=  0  for  z  56 

where 

n  =  INT(F')  for  i  =  1 
=  INT(F')  +  1  for  i  =  2. 

and  where  a  ,  b  , - ,  r  are  given  in  Table  Cl  for  n  =  0(1)5. 

n  n  n  s 

Then  findy(F')  =  rj  >  (t?.,  -  q  ^  ’  |  F'  -  INT(F')}. 

Find  P  corresponding  to  v(F'). 

The  recommended  formula  is  the  same  as  in  Appendix  B. 

For  0  f  F  ?  5,  set  PL(y},  s',  F,  r)~  P 

For  5  <  F  1  10,  set  PL  (y  ,  s’,  F,  r)  -  1  -  P. 


<C3) 
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The  following  BASIC  program  (filename  B:PYOSPER.  BAS)  estimates 
PL(Yq.  s'.  F,  r).  It  requires  these  inputs: 

seale  distance,  entered  as  R  (km), 

line  length,  entered  as  SP  (km), 

fraction  of  the  line,  entered  as  F  (tenths), 

the  cumulative  probability  for  the  threshold  value  (X),  entered  as  PO. 
and  yields  the  output  (P),  the  probability  of  the  event. 


10  REM  B: PYOSPFR. BAS  ON  GRINGORTEN  DISK  HO. 4.  IB  DEC  1986 

15  REM 

20  REM  THE  PROBABILITY  THAT  X.  OR  ITS  END  (YO!  WHOSE  CLIMATIC  FREQUENCY  IS  PO. 
WILL  BE  EXCEEDED  IN  (FROTHS  OR  LESS  OF  A  LINE  ( SP  KM):  PARAMETER  IS  SCALE 
DISTANCE  <R> 

30  REM  PO  HEREIN  WILL  BE  CUMULATIVE  PROBABILITY 

40  REM  THE  FOLLOWING  ARE  CONSTANTS  FOR  THE  ALGORITHMS 

50  AO  =  2.30753 

60  A1  =  . 27061 

70  B1  =  .99229 

80  B2  =  .04481 

90  Cl  *  .196854 

100  C2  *  .115194 

110  C3  =  .000344 

120  C4  =  . 019527 

130  DIM  AN(6).EN(6).CN(6) 

131  DIM  DN ( 6 ) .  FSN ( 6 ) .  GN<6> 

132  DIM  HN ( 6 ) . ON ( 6 ) .  RN  (  6 ) 

180  FOR  I  =  0  TO  5 

190  READ  AN( I ) , BN < I > . CN < I ) 

191  READ  DN( I ) . FSN( I  > . GN( I > 

192  READ  HN( I ) , 0N( I ) . HN ( I  ) 

200  NEXT  I 

210  DATA  . 03, -.0366, .51. 1.0049. 7. 44E-4.  . 82. 6. 06E-3. -1. 44E-5. 1.  02 

220  DATA  . 0346. -. 0366, . 47.  1.  0034.  5.  08E-4, .87.1. 89E-3.  -3. T^E - 7. 1. 46 

230  DATA  . 0272, -. 0283,  .  46.  1.0054.  2.  7~E-4. . 95. 2. 69E-3,  -2. 4SE-8.  1.  ei 

240  DATA  . 0207. -. 0213.  .  44.  1.  0058.  2.  74E-4. . 95.  -9.  76E-4.  -2. 34E-12. 2. 93 

250  DATA  . 0096, -. 0103.  .  44.  1. 0021.  3. ?2E-4.  . 91.  -2.  04E-3.  -7. 21E-17. 4.  17 

260  DATA  .  0000.. 0000.  . 00 . 1 . 0022. 3. 53E - 4 . . 92 . . 0000 . . 0 000 . . 00 

340  REM  STEP  0  FOR  THE  INPUTS 

350  INPUT  "SCALE  DISTANCE* “ :  R 

360  INPUT  "LINE  LENGTH* " :  SP 

370  INPUT  "FRACTION  OF  LINE  IN  TENTHS*": F 

380  INPUT;  "PO=": PO 

390  REM  STEP  1  TO  TRANSFORM  THE  VAR  I AELE  INTO  ALGORITHM  ENTRIES 

400  REM  TO  CALCULATE  2 

410  TEMP  =  SP/R 

420  Z  =  LOG ( TEMP ) /LOG  (  2  ) 

430  REM  TO  CALCULATE  YO  FROM  PO 
440  IF  POO.  5  THEN  480 
450  K  *  -1 

460  T  =  SQR ( LOG  < 1/PO'  2  )  ) 

470  GOTO  500 
480  K  =  1 

490  T  =  SQR ( LOG ( 1 '  ( 1 -PO  > ' 2 )  ) 

500  YO  =  K*(T-(A0  +  A1»T)/<1  +  B1«T*B2*T-  2M 

510  REM  TO  SUBSTITUTE  YOP  =  -YO  AND  FP  =  10 -F  WHEN  F',.5 

520  IF  F05  THEN  560 

530  FP  =  F 

540  YOP  =  YO 

550  GOTO  580 

560  FP  =  10-F 

570  YOP  =  -YO 

580  REM  STEP  2  TO  CALCULATE  THE  END  <Y>  OF  THE  DESIRED  FROEAEILITY 

589  REM  WHEN  FP=5  IT  SHOULD  EE  LESSENED  BY  A  TRIFLE.  TO  AVOID  A  ROUNDOFF  OF  ERRO 

R  IN  BASIC 

590  IF  FP  =  5  THEN  FP  =  4.999999 

600  FLOWER  *  INT(FP) 

610  FUPPE?  *  FLOWER  *■  1 

620  FOR  FI  =  FLOWER  TO  FUPPER 


83 


630  I  =  FI 

640  A  =  AN  < I  )  +  BN<I>«EXP<CN1  I)«Z> 

650  B  =  DN ( I )  *FSN( I >  »  EXP ( GN  < I ) *2) 

660  IF  2  =  < 6  THEN  700 
665  HEM  WHEN  Z>6 

670  C  =  HN ( I )  +  GN( I ) *EXP  <  SN ( I )  »Z> 

690  GOTO  710 
700  C=0 

710  IF  FI  =  FLOWER  THEN  K=1 
720  IF  FI  =  FUPPER  THEN  K  =  2 
730  ETA(K)  *  A  *■  E*YOP  ♦  C»Y0P'2 
740  NEXT  FI 

745  YP  =  ETA ( 1 )  +  (£TA(2)  -  ETA < 1 >)«< FP -FLOWER ) 

750  IF  F'-5  THEN  Y  =  -YP 
760  IF  F  =  05  THEN  Y  =  YP 

770  HEM  STEP  3  TO  CALCULATE  THE  DESIRED  PROBABILITY 

780  IF  Y>0  THEN  830 

790  L  «  0 

800  M  =  1 

810  X  =  -Y 

820  GOTO  870 

830  REM  WHEN  Y  ^0 

840  L=1 

850  M= -1 

860  X  =  Y 

870  DEN  =  (  1  +  C 1  »X +  C2* X'"2  +  C 3« X  3  +  C4 » X  4  )  '  4 
880  P  =  L  ♦  M / ( 2* DEN  > 

890  PRINT  "  THEN  P  =  “ ;  P 
900  GOTO  380 
1000  END 


Appendix  D 

Eitimatni  tin  Probability  That  a  ThrashaM  Valaa  (y0) 
Will  Nat  ba  Eiciatat  Ovar  A  ay  Liaa  of  Laaath  (s') 
Within  a  Longar  Liaa  of  Traval  (T) 


Purpose:  To  estimate  PI(yQ.  s',  T,  r),  the  probability  that  (yQ)  is  the  lowest 

maximum  END  in  a  linear  interval  (s'  km)  somewhere  along  a  line  of 
travel  of  length  (T  km);  scale-distance  parameter  is  r  km. 

Algorithm: 

Step  0.  Begin  with 

Single-point  cumulative  probability  (PQ)  of  a  threshold  condition 
(X),  with  corresponding  END  (yQ); 

Scale  distance  (r  km); 

Line  interval  (s'  km),  and; 

Overall  line  of  travel  (T  km). 

Step  1.  Compute  entries:  yQ,  z;  w,  as  needed. 

Find  z  =  ln(s' /  r)/ln2; 
u>  =  ln(T/s')/ln2. 

(Note:  Neither  z  nor  w  need  be  an  integer. ) 

Find  y  corresponding  to  Po.  The  recommended  formula  is  the  same  as  in 
Appendix  A. 

Step  2.  To  solve  for  the  probability: 

PI  (y  ,  z,  w)  =  P  (y) 
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find  y  (y  z,  w)  by  interpolation: 

y(yo,z,o>)  =  y  (yQ.  z?.  u>,).  (l-FRA(z)).  (l-FRA(w)) 

+  y  <yo.  Zp»  (l-FRA(z)).  (FRA(u>» 

+  y(yQ,  zu.  uJj).  (FRA(z)).  (l-FRA(u))) 

+  y(yo.  zu,  aO.  (FRA(Z)).(FRAM) 

where 

z?  =  INT  (z) 

zu  =  INT  (z)  +  1 

u >£  =  INT  M 

odu  =  INT  (a))  *  1 

FRA(z)  =  z  -  INT  (z) 

FRA(co)  =  a>  -  INT  M 

and  where 

y  <yo«  =  A  (i,  j)  +  B(i,  j)yo  4  C  (i,  i )  y 0 2 

where  A  (i,  j),  B(i,  j),  C  (i,  j)  are  given  in  Tables  Dl,  D2,  D3 
respectively. 

Lastly,  find  P  corresponding  to  y(y  ;  z,  u;).  The  recommended  formula  is  the  same 
as  in  Appendix  A. 

Set  PI  (yQ,  s',  T,  r)  =  P 
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Table  Dl.  A  (i,  j)  for  Given  z.,  ov. 


2  3  4  56  78  9  10  11 

0.025  0.  06  0.  13  0.  25  0.46  0.82  1.  35  2.  09  2.  95  3.62 

0.05  0.  12  0.  23  0.43  0.75  1.  26  2.  09  3.  09  4.50 

0.09  0.  20  0.40  0.72  1.  22  2.  02  3.  14  4.  79 

0.  14  0.  34  0.66  1.  18  1.  97  2.  95  4.52 

0.20  0.51  1.  03  1.84  2.84  4.  26 

0.29  0.  81  1.55  2.  52  3.  93 

0.37  1.11  2.01  3.30 
0.39  1.  23  2.  29 

0.13  0.94 

-0.78 


0.  125  -0,01 
05 


Table  D2.  B(i,  j)  for  Given  z 


Table  D3.  C  (i,  j)  for  Given  z^,  ac 


♦ 


I 


The  following  BASIC  program  (filename  B:PYOSPTR.  BAS)  estimates  PI  (y  ,  s’,  T,  r). 
It  requires  these  inputs: 

(a)  the  single-point  probability  of  exceeding  the  threshold  (X),  entered 
as  POP.  It  is  changed  by  a  step  in  the  program  to  cumulative 
probability  (PO). 

(b)  scale  distance,  entered  as  R  (km), 

(c)  the  overall  length  of  the  line  of  travel,  entered  as  T  (km), 

(d)  the  length  of  the  window,  entered  as  SP  (km), 

and  yields  the  output  (p),  the  probability  of  the  event. 
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10  REM  B: PYOSPTR. EAS  ON  GRINGORTEN  DISK  NO. 5.  10  MAR  1*87 

11  PRINT 

20  REM  BASIC  PROGRAM  FOR  APPENDIX  D  OF  REPORT  ON  3D-ESU  MODEL 
'30  REM  TO  FIND  THE  PROBABILITY  OF  THE  LOWEST  MAXIMUM  <  Y  0  >  IN  AN  INTERVAL  (S'! 
WITHIN  A  LARGER  INTERVALS) 

4  0  -  D I  M  A  (  1  0  .  1 4  ) 

50  DIM  B ( 1 0 . 1 4 ) 

60  DIM  0(10.14) 

70  FOR  I  =  0  TO  9 
80  2  =  I  -  1 
85  TWMINZ  =  12-Z 
90  FOR  J  =  0  TO  TWMINZ 


1 0  0 

READ 

At  I. J) 

110 

NEXT 

J 

120 

NEXT 

I 

121 

DATA 

. 01, . 02, . 0 

25, . 06, . 13.  .  25,  . 4b,  . 

82,  1. 35, 2.  09. 2. 95. 3. 62,5. 32 

122 

DATA 

- .  0 1 .  .  0 1 .  . 

05, . 12, . 23,  .  43.  .75, 1 

.  26,  2.  09.3. 09, 4. 5, 6. 26, 6. 87 

123 

DATA 

-.02,  .  02,  . 

09,  . 2, . 4, . 72,  1. 22.  2. 

02. 3. 14. 4.79.6. 62.  8.6 

124 

DATA 

-.05, . 02, . 

14,  .  34.  .  66,  1. 18, 1. 97 

,  2. 95, 4.52,6. 33, 8. 65 

125 

DATA 

-.  125, -. 01 

,  .  2,  .  51,  1.  03.  1.  84,  2. 

84, 4. 26, 6. 3,7. 74 

126 

DATA 

-.  25.  -. 05, 

. 29,  .81,1.55. 2. 52. 3. 

93. 5. 82, 8. 18 

127 

DATA 

-. 45, -. 14, 

.  37,  1.  11,  2.  01,  3.  3,  5. 

03,  7.91 

128 

DATA 

-.76,  -.31, 

. 39,  1.  ?3,  2.  29, 3. 45, 5 

.  35 

129 

DATA 

-1. 23, -.62 

, . 13, . 94. 1 . 9. 3. 16 

130 

DATA 

-1. 91.  -1.  2 

8. -. 78. -. 097, . 54 

131 

FOR  I 

=  0  TO  * 

140 

Z  =  I 

-1 

145 

THMIN 

Z  =  12  -  Z 

150 

FOR  J 

=  0  TO  TWMINZ 

.. 

160 

READ 

E(I. J) 

170 

NEXT 

J 

180 

NEXT 

I 

190 

DATA 

1.  1,  1,  1.  1.  01.  1.  03. 

1.06.1.12.1.24.1.44. 

1. 

65.  1. 6B, 2. 51. 3. 3 

191 

DATA 

1,  1.  1,  1.  01,  1.  02,  1. 

06, 1.11,1. 26, 1. 54, 1. 

85 

, 2. 44. 3. 13, 3. 04 

192 

DATA 

. 99. 1, 1, 1. 02, 1. 06, 

1.  12,  1.  24,  1. 5,  1. 96,  2 

.  f 

2,3.  41,4.  14 

193 

DATA 

1,1.01.1.01,1. 06.  1 

.  10,  1. 2,  1. 51, 1. 84, 2. 

55 

, 3. 34. 4. 35 

194 

DATA 

1,  .  99,  1.  04.  1.  11,  1. 

23, 1. 5,  1.  83,  2. 39, 3. 3 

1, 

3.61 

195 

DATA 

1.  01, 1. 02, 1. 1,  1.  23 

,  1. 42, 1. 8,  2.  42. 3. 31. 

4. 

38 

196 

DATA 

1.  03, 1. 07, 1. 2, 1. 44 

,  1. 77,. 2.  39.  3.  27.  4.97 

197 

DATA 

1.  1,  1.  15,  1. 39,  1.  71 

,  2. 23, 2. 71. 3.91 

198 

DATA 

1.  21,  1. 33,  1. 63,  2.  2 

,  2. 87, 3. 99 

199 

DATA 

1. 45,  1. 65,  2, 2. 52. 3. 6 

200 

FOR  I 

=  0  TO  9 

210 

Z  =  I 

-  1 

220 

TWMIN 

Z  =  12-Z 

230 

FOR  J 

=  0  TO  TUMINZ 

24  0 

READ 

C  (  I ,  J  ) 

250 

NEXT 

J 

26  0 

NEXT 

I 

270 

DATA 

0 ,  0 , 0  ,  0 , 0  ,  . 

0 1 ,  .  0 1 

,  0 ,  .  0 1 .  .  0  2  ,  .  0  3 ,  0 .  . 

12. . 22 

271 

DATA 

0, 0,  0,  0,  0,  . 

0 1 ,  .  0 1 

,.02,. 05,. 08,. 15,. 

23, . 1? 

2  72 

DATA 

0 ,  0 ,  0 ,  .  0 1 ,  . 

01, .01 

,  .  0  2  ,  .  0  4 ,  .  1 .  .  2 1 ,  .  2 

7.  .  35 

273 

DATA 

o 

1^1 

o 

>-* 

01.  0,  . 

05, . 08 ,  .  17, . 27,  .  39 

274 

DATA 

0,  . 01,  . 01. . 

02, . 02 

,  .  04,  .  07,  .  13,  . 25.  . 

?  c, 

275 

DATA 

.  01.  . 01. . 01 

,  .02,  . 

03, . 07, . 15, . 27, . 4 

276 

DATA 

.  0 1 .  .  0 1 ,  .  0  2 

,  .  03,  . 

07, . 16,  . 29, . 57 

277 

DATA 

0 ,  0  ,  .  0  3 ,  ,  0  5 

r  -14,  . 

21, . 45 

270 

DATA 

0 ,  0 ,  .  0  5 ,  .  1  0 

,  .32.  . 

6 

279 

DATA 

-  .  0  3 ,  - .  0  8 ,  . 

04, . 14 

,  .75 

300 

INPUT 

1  "PR0EAEILITY  OF 

CLOUD  COVER”:  POP 

% 


4 


V 


00 


< 


■\  * 


301  REM  PROBABILITY  OF  CFLOS  IS  PO  WHERE 
310  PO  =  1  -  POP 

320  REM  TO  FIND  YO.  THE  END  OF  PO 
330  GOSUE  700 

340  INPUT  “SCALE  DISTANCE  =  R" :  R 
350  INPUT  “OVERALL  DISTANCED"  sT 

359  LIMIT  =  256*R:  IF  TCLIMIT  THEN  LIMIT  =  T 

360  PRINT  “LIMIT  WINDOW  (SP)  TO";  <T/2  13i:  "  :  SP  :  “ : LIM IT 

361  INPUT: “WHEN  SP=”: SP 
370  S  =  SP/R 

300  Z-  LOG ( S ) /LOG ( 2 ) 

301  IF  22-1  THEN  304 

302  REM  WHEN  2  =  2-1 

303  Z  =  -l 

304  REM  WHEN  Z>-1 

305  IF  2<B  THEN  390 

306  REM  WHEN  228 

387  PRINT  “  NO  ANSWER:  ALGORITHM  WAS  NOT  EXTENDED  FAR  ENOUGH. 
380  GOTO  361 

390  OMEGA  *  LOG < T/ SP ) /LOG < 2 > 

400  ZL  =  INT ( Z ) 

410  ZU  =  INT ( Z)  +1 
420  ML  *  INT ( OMEGA ) 

430  HU  *  ML  +  1 

435  REM  FOR  THE  FIRST  TERM  (Yl) 

440  ZP  =  ZL 
450  HP  =  ML 

460  GOSUB  810 

461  FRAZ  =  Z  -  INT (Z) 

462  FRAH  =  ONEGA  -  ML 

470  Yl  *  YP* ( 1-FRAZ  >  * ( 1-FRAW  > 

480  REM  FOR  THE  SECOND  TERM  (Y2) 

4  )<>  ZP  *  ZL 
500  HP  *  HU 
510  GOSUB  810 

520  Y2  =  YP* ( 1-FRAZ )  t(FRAH) 

530  REM  FOR  THE  THIRD  TERM  (Y3> 

540  ZP  =  ZU 
550  HP  -  HL 
560  GOSUB  810 

570  Y3  *  YP* ( FRAZ ) « ( 1-FRAH ) 

600  REM  FOR  THE  4TH  TERM  (Y4) 

610  ZP  *  ZU 

620  HP  *  HU 

630  GOSUB  810 

640  Y4  *  YP*  < FRAZ ) « ( FRAH ) 

650  Y  *  YI+Y2+Y3+Y4 
660  GOSUE  860 

670  PRINT  ’  THE  PROBABILITY  =  “:P 
675  GOTO  361 

700  REM  THE  SUBROUTINE  TO  FIND  YO  GIVEN  PO 

701  AO  «  2.30753 

702  A1  ■  . 27061 

703  B1  *  .99229 

704  B2  >  .04481 

710  IF  P02.5  THEN  760 
720  REN  HHEN  PO  ><.5 
730  TP  ■  SOR (LOGt 1/P0A2 ) ) 

740  K  «-l 
750  GOTO  790 
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7*0  HEM  WHEN  P0>.5 

770  TP  =  SQR ( LOG  < 1 / ( 1 - P  0  )  '  2  )  ) 

780  K  =  1 

790  YO  =  K»<TP-<A0-*A1«TP)/<1*B1«TP*E2«TP'Z)> 

800  RETURN 

810  REM  SUBROUTINE  TO  FIND  THE  END  <Yi  FOR  GIVEN  Z.  OMEGA  AND  YC 
820  I  =  ZP  +  1 
830  J-  HP 

840  YP  =  A(I.J)  +  E  (  I ,  J  )  »  Y  0  +  CtIfJ)»Y0'2 

850  RETURN 

8*0  REM  SUBROUTINE  TO  FIND  THE  PROBABILITY  P  CORRESPONDING  TO  Y 

870  Cl  =  .19*854 

880  C  2  =  .115194 

890  C3  =  .000344 

900  C 4  =  .019527 

910  X  =  ABS(Y) 

915  DEN  =(1  ♦  C 1  * X * C 2 * X  2  ♦  C3*X'3  ♦  C4*X'4>  4 

920  IF  Y<0  THEN  980 

930  REM  HHEN  Y*>0 

950  L  =  1 

9*0  H  =-l 

970  GOTO  1010 

980  REM  HHEN  Y  20 

990  L  =  0 

1000  M  =  1 

1010  P  =  L  ♦  MM. 5/DEN) 

1015  RETURN 
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